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LATERAL DYNAMIC STIFFNESS OF PILES ar 


_ KEY WORDS: Dynamic response; Dynamics; Foundations; Lateral loads __ 


Giles); Piles; Soil | dynamics; Soil properties; Soil-structure 


: Stiffness; Vibration 


ABSTRACT: The effect of nonlinear ae nari! (without slippage of gapping) on the 
dynamic response of piles subjected to lateral loads is explored, using a finite element 
_ model for the soil region adjoining the pile, a consistent boundary matrix at some 
distance to reproduce radiation effects, and an iterative, 


~ technique to estimate the variation of soil properties with level of strain. Results 
- obtained with this procedure are compared with those that would result from 
application of the p-y curves and the differences are discussed. ig), pede 
_ REFERENCE: Angelides, Demosthenes C. (Sr. Consulting Struct. Engr., Research and 

_ Development Div., McDermott Inc., P.O. Box 60035, New Orleans, La. 70160), and 
- Roesset Jose M., “Nonlinear Lateral Dynamic Stiffness of Piles,” Journal of the 
“November Engineering Division, ASCE, Vol. 107, No. GT11, Proc. Paper 16635, 


16634 COUPLED AND WATER FLOWS BURIED 
KEY WORDS: Backfills; ‘Buried | cables; Computer analysis; Diffusion; Field fort 
| tests; Heat transfer; Numerical analysis; ‘Soil moisture Temperature 


distribution; Temperature gradients pa 
ABSTRACT: Theoretical and experimental tution: of heat flow and moisture 


_ movement under thermal gradients are described with particular emphasis on flows — 
_ around buried electrical transmission cables. A one-dimensional finite difference — 
computer program for analysis of coupled heat and water flows in layered systems 
based on the Philip and DeVries theory is developed. Theoretical analysis indicate that 
the most important parameters affecting rate and amount of moisture migration away | 
from buried heat sources are the sources temperature and the soil type. A field test — 
was conducted to evaluate the heat and moisture flow predictive methods. The results 
_ showed that water content predictions using the coupled flow computer program match — 
poorly with the actual values, and that predicted values are very sensitive to the soil | 
properties used. A finite element two-dimensional computer program is used to predict 
the variation of temperature distributions around the buried cable with time after a4 
initiation of cable heating. The results show very good agreement between measured 
REFERENCE: Abdel-Hadi, Omar N. (Geotechnical Engr. Geotechnical Consultants, 
-Inc., San Francisco, Calif.), and Mitchell, James K., “Coupled Heat and Water Flows 
- Around Buried Cables,” Journal of the Geotechnical Engineering Division, ASCE, Vol. 
107 No. GT11, Proc. Paper 16634, November, 1981, PP. 1461-1487 erg dt 9D 


_ ABSTRACT: Curent. practices and concepts for establishing the tolerable settlement 4 p 
of buildings are reviewed. The role of differential settlement is emphasized. Most 
current tolerance limits are shown to be based on the works of Skempton, MacDonald, eh 
Polshin, and Tokar. These criteria are compared with field observations of settlement — 
and damage. A simple beam analogy, which was proposed by Burland and Wroth, is — 
used to illustrate factors that influence tolerable settlement. The effects of the critical — 
tensile strain of the structural materials, the L/H ratio of the building, the relative \ 
stiffness of the structure in shear and bending, the longitudinal stiffness of the of 


foundation, and the shape of the settlement pattern are demonstrated. Different criteria — 
are required for different types of structures. The allowable eam is smallest for 
cases in which the settlement curve is concave downward. 


REFERENCE: Wahls, ‘Harvey E. (Prof. of Civ, Engrg. , North Carolina State Univ., g 
Raleigh, N. C.), “Tolerable Settlement of Buildings,” Journal of the Geotechnical — : 
Engineering Division, ASCE, Vol. 107, No. Proc. e's 16628, 
1981, pp. 1489-15000 
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16636 INSTABILITY OF AMUAY CLIFFSIDE HAS RE 


KEY WORDS: Case Clays; Cliffs; Field investigations; 


yt period indicate that rising pore pressures trigger slides along the cliff. Visual — 
> observation, surveys, and field exploration show that most of the landslides have 
_ wedge-type failure geometry. Shear occurs through a nearly horizontal stratum of weak 
clay, brown fat clay. The strength of the fat clay decreases with deformation. Fluids 
Beng ox by domestic and industrial activities have resulted in a perched water ofl 
a above the stratum of fat clay. For more than two decades, the pore pressure on top of 
the fat clay hes risen. A wedge type of stability analysis can give the correct factor of — 
safety — i.e., 1.0 — for the slides which have occurred. 
| REFERENCE: Lambe, T. William (Prof., Dept. of Civ. Engrg., , Massachusetts Inst. of — a 
Tech., Cambridge, Mass. 02139), Silva, Francisco, and Marr, WwW. Allen, “Instability of aw 
| Amuay Cliffside,” Journal of the Geotechnical Engineering Division, ASCE, Vet. 107, 4 
No. GT11, Proc. Paper 16636, November, 1981, pp. 1505- 15-1520 —_ 
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16663 FRICTION CAPACITY OF PILES DRIVEN INTO CLAY One. 
GAD 

-KEY WORDS: Clay | soils; Friction resistance; Penetration; Pile foundations; gO" we = 

Pile friction; Pile loading tests; ‘Static pile wae 

_ ABSTRACT: — studies on axial pile capacity in clays have shown that the 

average frictional resistance, expressed as a fraction of the average undrained shear 

_ Strength or average effective overburden pressure, decreases with increasing pile a5 a 

_ penetration. Procedures to compute shaft friction are reviewed, and the effect of pile 

«length on the development of shaft friction on piles in clay is examined in terms of the if 

relative pile-soil stiffness and lateral pile movements during installation. Correlations — 4 

| are developed to relate shaft friction coefficients a B A to pile length, relative pile-soil _ 

eS stiffness, and soil stress history. Procedures are recommended to compute the friction i 

+ REFERENCE: Kraft, Leland M. (Mgr., Special Projects Group, McClelland Engrs, a 

| Inc., Houston, Tex.), Focht, John A., Jr., and Amerasinghe, Srinath ———— t 

g Capacity of Piles Driven into Clay,” Journal of the Geotechnical Engineering Division, 
_ASCE, Vol. 107, No. GT11, Proc. Paper 16663, November, 1981, pp. 1521- MAbs tne 
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16653 THEORETICAL T-Z 30 07 we 


se settlement; Predicting Settlement Stress; Subgrades 

‘Most procedures to compute “the: ‘(stress- displacement) response for 
axial pile loading are empirical, based mostly on data for pile diameters less than 18 in. 
(0.15 m) and pile penetrations less than 100 ft (30 m). These procedures may not be 
appropriate for pile and soil conditions that differ from those on which the procedures — 
were developed. Therefore, a theoretical procedure is developed to generate t-z curves | 
for predicting pile movements under axial load. The procedure uses an a 
one-dimensional, elastic continuum approach to define the t-z response up to the 
maximum t response and laboratory simulation to define the post-peak response. The 
procedure was used with several case studies. Predicted response of the pile under 
axial loading was in good agreement with the measured response. = 
REFERENCE: Kraft, Leland M., Jr. (Mgr., Special Projects Group, i 
Engrs., Inc., Houston, Tex. 77081), Ray, Richard P., and Kagawa, Takaaki, 
“Theoretical t-z Curves,” Journal of the Geotechnical Engineering Division, _ ASCE, 
Vol. 107, No. GT11, Proc. heal 16653, November, 1981, pp. 1543-1561 
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SEEPAGE FROM FREEWATER A Oo 
KEY WORDS: Complex variables; Dams (earth); ne wastes; in 
Models; Seepage; Soil mechanics; Tailings 
ABSTRACT: A proposed model for seepage discharge into a porous etentuses from .Y 
free water above impermeable tailings is described. It is assumed that the discharge is r 4 
controlled by saturated seepage in an entry zone adjacent to the free water interface. - 
The isotropic entry zone has been analysed by complex variable methods and a general — 
solution obtained for seepage discharge quantities. Application of the solution to © 
anisotropic embankments is described. The significance of the theoretical results for — 
- design is discussed and a simple formula proposed for estimating the Fr, ~ ons 
ence. Isaacs, Lewis T. (Sr. Lect., Dept. of Civ. Engrg., Univ. of Queensland, . 
Brisbane, Australia), and Hunt, Bruce, “Seepage From Freewater Above —— 
Journal of the Geotechnical 1981 pp. 156 Division, ASCE, Vol. 107, No. 
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use has been published by the American Society for Testing and Materials. Copies of this — 


publication (ASTM E-380) can be puscheses | from ASCE at a price of $3.00 each; ime must — 


All authors of Journal papers are being asked to prepare their papers in this dual-unit format. 


To provide preliminary assistance to authors, the following list of tinned factors and oats 
recommended by the ASCE Committee on Metrication 
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meters (m) 

_ miles (miles) Woy kilometers km), 
Square inches (sq in. Square millimeters (mm’°) 
square feet(sqft) square meters (m 
"Square yards (sq square meters (m’) 
square miles (sq miles) square kilometers (km & 

cubic inches (cu in) cubic millimeters (mm 
“cubic feet (cu ft) cubic meters (m ~ 

pounds (Ib) mass kilograms (kg) 
(ton) mass kilograms (kg) 


newtons (N) 


newtons (N) 
pounds per square (psf) pascals (Pa) 
pounds per square inch (psi) (kPa) 89 
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By! Demosthenes C. A. M. and ASCE 
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asviu9 adi lo .mot Heb geod) dtiw 
everal procedures have been developed in recent years to determine the ~ 

7 - dynamic response of piles subjected to horizontal or vertical loads of relatively — 
small amplitude, , and where the soil can be assumed to behave a a ec 


elastic material. Most of these theories yield very similar results eines that 
the key features of the problem are now | well understood. ae 


soil behavior, slippage, and eventual gapping, play a fundamental role in > 
response of piles to cyclic loads of moderate to Lye amplitude. These charac- 
teristics ha ve been implicitly incorporated i in models of piles as beam on elas 
foundation where the Winkler type springs are given by the p -y or f-z curves. 
_ The concept of the p-y curves was introduced initially by McClelland and a 


Focht (12) as the first significant pitames to model nonlinear behavior in the 
his approach using a finite ‘difference formulation | to discretize ‘the pile and 
the nonlinear foundation. Matlock (11) presented p-y curves and their analytical 
expressions for soft clays, and Reese, et. al., (15,16) established similar curves © 
for stiff clays. A distinction was made by Matlock (11) between static (implying 
monotonic) and cyclic loading. The cyclic loading p-y curves deviate from the 
static ones when the deflection at the mudline reaches a limiting value due 
to cyclic degradation of the soil. It should be noted that these cyclic p-y curves © 
represent an envelope of the cyclic behavior, and do not describe explicitly 
the reduction of stiffness of the soil as a function of the ‘magnitude of the 


force for a transient excitation, nor as a function of the number of oS 


More recently, Stevens and Audibert (18), by analyzing previous experimental 4 


"information from several sources, § Suggested | modification of the p-y curves given 
by Matlock (11). They: recommended a different relation between the characteristic ‘ 


a Consulting Struct. Engr., Research and Development Div., Inc., Py oO. 
 *Prof. of Civ. Engrg., The Univ. of Texas at Austin, Austin, Tex. 78712, NS ilo 
_ Note.—Discussion open until April 1, 1982. To extend the closing date one e month, 
a written request must be filed with the Manager of Technical and Professional Publications, 
i. ASCE. Manuscript was submitted for review for possible publication on February 25, 
1981. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings _ 
of the American Society of Civil Engineers, OASCE, Vol. 107, No. GT11, November, : 
1981. ISSN 0093-6405 /81 /0011-1443 /$01.00. 
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soil deflection, y. , the deflection at one half the ultimate soil resistance, 
— the pile pg ‘ead a modification of the ultimate lateral soil ‘sevetance 


_ Methods of ‘dynamic pile analysis based on p-y curves can tone ecctiunndiies 
“quite realistically; they cannot directly handle, however, soil inertia or radiation 

In this paper an attempt is made to explore the effect of nonlinear soil behavior 

(without slippage or gapping) on the dynamic response of piles: subjected _ 
lateral loads, using a finite element model for the soil region adjoining the i 
Fi pile, a consistent boundary matrix at some distance to reproduce radiation effects, 

and an interative equivalent linearization technique to estimate the variation 

of soil properties with level of strain. Results obtained with this procedure 

are compared with those derived from application « of the ‘py curves and the 


The main model used in this study is an extension | of the formulation presented __ 
by Blaney, | et. al, (4) for the linear elastic case. A cylindrical region of soil 
4 
surrounding the pile (Fig. 1) is discretized ‘using toroidal finite elements. The Ps 
consistent boundary matrix derived by Kausel (7) for problems in cylindrical 
_ coordinates is placed at the edge of the core region. The pile is modeled as - 
a series of beam segments attached through rigid links to the finite — 
An iterative equivalent | linearization | is ap plied to the ‘elements of the core 
region, following the same general method used by “Seed and Idriss’ (17) for 
_ Seismic studies of soil deposits (one-dimensional pure shear condition), and 
7 implemented in such well known programs as LUSH (8) and FLUSH (9) for 
two- dimensional or pseudo- -three-dimensional conditions. In each cycle of analy-— = 
sis, soil Properties: are assumed to be constant, with a shear modulus, G, , and Py 
a hysteretic damping | ratio, | a consistent with the amplitude of the shear strain, 
Yas obtained as a result of the previous cycle. The selection of this strain 
is clearly arbitrary. Moreover, the value of this strain changes i in the circumferen- gfe 


‘tial direction with the sine of the angle 6, but the model | assumed a = 


e set of values of G and D, for each toroidal element, _ making rage) ba 
independent of 8 and equal to the values that would occur at 6 =+ orl a 
- In addition to the finite element model, special purpose ‘computer programs 
- were developed discretizing the pile with linear members and replacing the soil 
: with lumped springs at the nodes. The spring characteristics were derived som) 
‘Matlock’ s and Stevens and Audibert’ s p-y curves. For both models, hollow, 
steel piles with aspect ratio (length, g* over | diameter, d) equal to 45, diameters | 
ft, 2 ft and ft, and thicknesses of 2. Sin., 1.25 i in., and 0. 625i in.., ‘fedpectively, 
were considered. These result in a value of K, = E,1/E, L* of 1.82 x 10° 
for all cases, with E, the pile modulus of elasticity E. , the soil modulus of | 
L the length of pile, and /its moment of inertia. 
Since the objective of this study was to observe trends i in behavior, no attempt 
pe made to simulate any particular soil profile. Even so, it was desirable 
: have soil properties reasonably typical of soft clays. One of the problems 
encountered in comparing results of analytical or numerical formulations based > 


| 
| 

i 

| 
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LATERAL AL DYNAMIC STIFFNESS 
nose of p p-y curves is to select consistent wae 
_ soil parameters. For the finite element solution, the soil behavior is charac- 
an terized by the shear modulus, and its variation with shear strain; a second | 
_ elastic constant, such as the bulk modulus and the total unit weight, is required | 
_ compute the shear wave velocity. The p-y curves, on the other hand, are 
normally developed from the undrained shear r strength, the effective. (submerged) 
- as weight of the soil, and a a characteristics Strain, €,,, which is equal to the 
strain at one-half the maximum deviator stress in an undrained laboratory triaxial 
test. The relationship between shear modulus and undrained shear strength is 
easily established. ‘Tt was decided, therefore, to use a soil for | which 
8,135 EN fmm’ oF 28,477 pet}, 9.6 (96.14 kN 
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FIG. 2.—Force- Seteeten Relations along Pile—Comparison with Matlock, Matlock, 1970 (1 ow 
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FIG. . 3.—Force-Deflection Relations with Stevens and Audi- 


bert, 1979 (1 Ibf/ft = 14.59 N/m; 1 ft = 0.305 m 
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LATERAL OVNAGC C STIFFNESS 


experimental was The desired were the undrained 
shear strength, the dynamic shear modulus at small strains, and the variation 
of this modulus with shear strain amplitude, number of cycles, and frequency — 

(or period) of the excitation. The soil selected is a plastic, overconsolidated _ 
clay with an overconsolidation ratio of 4, extensively studied by NGI (13) and — 
q “4 Andersen (1,2). This experimental information provided the variation of shear 


ge with number of cycles, N, , and shear strain amplitude, y. To account 5 


for the effect of frequency, data reported by Fisher, et. al, (5) for a plastic 
ee clay with the same overconsolidation ratio and very similar characteristics were 
‘The soil properties used for the analysis were: modulus Gar, = 1,600 t/ m? ¢ 
, (5, 734 kN /m’ or 328,477 psf); undrained shear strength 9.8 t/m7 (96.14 kN/m? 
or: 2,007.06 psf); Poisson’ s ratio v = 0.49; variation of modulus G with shear 
Son Se eat) OH 
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d=2ft, t= 1.25in. 
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FIG. 4. —Real Part of Horizontal Soil aceon Versus Depth for Different Levels. of 
Force F (kips) [1 (Ibf/ft)/ft = 47.84 (N/m)/m; 1 0.305 
‘Strain y, number of cycles, N, and period of excitation, re” 


+ ay 
a given 
in which a = 216.16; b = 0.1038; ¢ = 2.777988; d = — : ates rotenges log = 
‘then natural logarithm. (The a above relation was obtained by fitting the experimental — ; 
curves given in Refs. 13, 2, 1, and 5. Details on n the achieved accuracy in 
these fittings are given in Ref. 3.) It was assumed that these properties were 
uniform over the soil deposit, 300 ft deep. The value of y used in the finite 
q element program is that of y,, previously explained. 
_ Following Hardin and Drnevich (6), the relation be between the - hysteretic damping ¥ 
ratio, D,, and the shear modulus was taken as 


d=4ft ,t= 2.5in. ,f=0.5 Hz ,N=lOcycles 
od) Yo 


solution 


008, 


t= 25in, 


- F.E.M. solution 
+i 


1 


gai with p-y cyclic 
200 400 600 800 
HORIZONTAL FORCE (kips) 


of with Level of Force, for Different Frequencies 
14.59 N/m; 1 kip = 4.45 KN; 1 ft = 0.305 m) eelubom vesde ort 


tmz 200 400 600 800 00 
tm: HORIZONTAL FORCE. (kips) 


No information was available on the value of Dus _ for the soil selected. =a 
4 Hardin and Drnevich suggest for saturated cohesive “soils ; an expression for 


4 D.,ax Which is a function of the effective mean principal stress, Tos the —: 


of cycles, N, and the frequency of ad 


A pe 1/2 
Pax = 031 + 0.015/ 0.0065 log N. 


any 
in 1 which = kg/c Since the modulus hed constant 
~ over the depth of the soil layer, it was decided to take also a constant friveriverine— 
the effect of for the 


Ja? ‘ 
In the finite element solution, forces” and or in the soil along, 


4ft t= ,N=10 cycles 


Fehorizontal force at top of pile pile (hips) 


are ate a discrete number of To establish a a comparison 
with the p-y curves, the stiffness matrix x of the soil cavity, ‘corresponding to 
a given level of forces at the head of the pile (computed with the final iterated 


- _ properties), was multiplied by the nodal displacements along the pile to 


obtain equivalent soil forces. The resulting force deformation relations are then — 
a function of the deformed shape and, therefore, of the relative stiffness of 


the pile with resp respect tothe soil. @ fecwr of When 
a. Fig. 2, the force-deflection relations obtained with this - procedure are 


q 
i 
Pow 


presented and compared to the p-y curves, derived from Matlock’s (11) rules, © 4 


4 for three pile diameters and three locations along the piles. The forces secehting 
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even larger when degradation starts to take place. The limiting 
a at which the cyclic p py curves deviate from 


of forces as the pile decreases. 
Several reasons may cause the difference in results: = 

ta . The iterative procedure is, by itself, only an approximation and an even 4 
- 2 one when dealing, as in the present case, with a truly three-dimensional : 
y _ 2. The assumed constitutive equations for the soil are only approximate, and 
the relation between stiffness (shear modules) and Strength is an additional 


Hz 


to snot iotnosiiot =F | | 


=4ft,t= ,t=2.5in. ,N=lO cycles 


400 600 800 1000 


HORIZONTAL FORCE (kips) 


, with Level of Force, for Different Frequencies (1 ibf/ft = 
14.59 kip = 4.45 kN; 1 ft = 0.305 


3. Results are based on 10 cycles of loading, but the p-y curves are supposed 
te represent a lower bound to the response after an unspecified number of 
a 4. Slippage and flow of soil immediately. adjacent to the f pile are not considered. = 
Studies with two-dimensional models by Yegian and Wright (20), including these 
effects, showed better agreement with the p-y curves, but still not a perfect | 
qa: The f fact that the soil cannot take tension near the surface and that gaps 
nd will open is also ignored. Studies by Thompson (19), with a _ two-dimensional 
me showed that the soil forces could decrease by a factor of two a 
= Fig. 3 shows, on the other hand, a comparison of the force-deflection relations — 
for the soil, obtained from the finite Gentemt model, with th the p- sie suggested 


— 


= 


by Stevens ‘Audibert (18). The agreement is now much better near the 
_ surface, but the forces corresponding to a given deflection from the p-y curves — 
are larger than those predicted by the finite element solution at some depth. P 
_ The finite element results seem to be between the p-y curves of Matlock and 

. _ The variation of the real part of the soil springs, k,, with depth, is shown 

in Fig. 4 for different levels of horizontal force at the pile head and for pile 
diameters of 4 ft, 2 ft, and 1 ft. It should be noted a. while there is nonlinear 


\ 


linear 
linear, 


d=4ft,t=2.5in.,N=10 cycles 
F= horizontal force of pile 


: FIG. 9. —Variation of c,, with Frequency, for Different Levels of Force (1 Ibf/ft = — 
14.59 N/m; 1 kip = 
f) Hil! 1? to isved fiw o> te 
figure. It appears, thus, that some t basic behavi 
by the nonlinear model. 
_ The finite element model, and the programs a using the p-y curves, 
were then used to investigate trends in horizontal pile stiffness and damping 
_ characteristics. The lateral stiffness of the pile was obtained in both cases by 
_ applying a force at the pile head, which was free to rotate, computing the © 
deflection through a series of iterations, and dividing the force by the top 
- deflection. Clearly this represents only one particular form of excitation (horizon- | 
tal force at the pile top, no moment or axial load). Some analyses were performed 
for more general conditions, such as fixed head piles with moments and axial _ 
forces, in addition to the lateral load; but no parametric studies were conducted py , 
for these cases, due to the large number of parameters involved. Results presented 
= are therefore limited to the case of a pile hinged at the top. dn, ten wt 6 


, is particularly large over a small depth. This depth seems to increase slightly - ; 
with level of force, and is comparable to the depth of reduced resistance, 
j 
4 
OS 
| 
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a? vary with frequency. The imaginary part, c,,, represents the loss of energy. * 
The ratio c,, / 2k. can then be considered as the effective damping 4 
at each fre requency. The D value independent « of frequency represents a hysteretic " 
: type damping associated with nonlinear soil behavior; a D value, increasing 
with frequency, represents a viscous type damping characteristic of 


The variations of k,, with level of horizontal force at the pile head, predicted 
- by the finite element solution, and d by the model with lumped nonlinear r springs © 7 


= on Matlock’ wey: ‘are shown in Fig. 5 for a pile diameter of 


< 
d= 4ft,t= .2.5in.,N=lOcycles 


“Ae 


Pat. 
OLS 1 200 600 800 1000 


4 ft and three values of the ¢ excitation frequency (0.5, 2, and 3 Hz). — curves — 
ar 
q with increasing | force. It should be noted, i however, that the level of force 
to initiate cyclic degradation with the p-y curves decreases with i increasing 
frequency. The stiffnesses predicted by the finite element solution are 2-2.5_ 
J times larger than those obtained with the p-y model, a result consistent with — 
that obtained for the equivalent soil springs. 
ott The variation of pile stiffness, k,,, with frequency for various force amplitudes, 
is further shown in Fig. 6 for the same pile diameter and the finite element _ 
- model. The slope of these curves is essentially invariant with the level of force, 
but the relative decrease of stiffness with frequency is larger as the force at 


| 
j 


‘The effect of the one ‘diameter is shown in The variation the 
as stiffness, with level of force, is shown for three pile diameters (4 ft, 2 
ft, and | ft), a frequency of 1 Hz and the two models (finite elements and 
lumped springs). Comparison of the results from the two formulations lends an 
the same observations made in relation to Fig. the stiffness predicted 
: by the finite element solution is 2-2.5 times larger than that obtained from 
_ The variation of the imaginary part, c,,, of the lateral pile stiffness, ~~ 
_ level of force at the pile head, is shown in Fig. 8 for a pile diameter of 4 
ft and frequencies of 0.5 Hz, 2 Hz, and 3 Hz. For the smallest frequency 


incor, 


| 
4 F=horizontal ‘force at top of pile (kips) 


- 1b 5 20 25 30 

FIG. 11.—Variation of c,, /2k,, with Frequency, for Different Levels bat Force (1 kip 


= 4.45 kN; 1 ft = 
45 KN; 1 ft = 0.305 m) ~of veto? Ys level te 


Hz), the c,, term increases initially with increasing force, reaches a maximum, 


then decreases. For the higher frequencies, it decreases monotonically. This = 
behavior can be explained by considering that there is an increase in internal a 
hysteretic type damping, due to soil nonlinearity, and a decrease in radiation 
damping, due to the softening of the core region with respect to the far field. 
For low frequencies, the increase in internal damping is larger than the decrease 
in radiation, up to a certain level of force. Beyond this level of force, the © 
internal damping remains essentially constant or increases very slowly, while 
the radiation damping keeps decreasing. For larger frequencies, the c,, term, 
due to radiation, becomes much more important than the contribution of the 

, decreases monotonically with 


7 
— 
4 
— increasing level OF force. These results correspond to the finite clement 
a . where the pile was assumed to remain in contact with the surrounding soil 
a = at all times. Due to the appearance of gaps, one might expect a reduction be 


The variation on of ‘with frequency, is s shown in n Fig. 9 for various levels 
_ of horizontal force at the pile head. For low frequencies, the lines cross each — 
* other. In fact, for frequencies smaller than the natural frequency of the soil | 
_ deposit (about 0.25 Hz in this case), there would be no radiation, and the only 
3 term w would be due to internal soil ques Over r the range of pape gen ie 


lOcycles 


10 


N=l0 cycles. 


800 
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FIG. 12.—Variation of c,, with Level of Force, for Different Pile Diameters rs (1 Ibf jt 
‘studied, term varies almost linearly, with the slope 
level of force i increases. of red! 
The above results are of interest if one intends to replace the pile by a - 
' spring with a stiffness, k,., and a viscous dashpot, c,, /w, both of which are . 
functions of the frequency of excitation, w(w = 2nf), and of the amplitude 
of of ‘the applied force force. It is a also interesting to consider the effective damping 
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ratio, taking i into account the variation of both terms. Fig. 10 sh 
m the variation of this term with level of force for the 4-ft diam pile, and the 
or three frequencies 0.5 Hz, 2 Hz, and 3 Hz. In all cases, the ratio mrs 


with increasing level of force. This indicates that the ‘reduction of the real 
tiffness, k,,, more than compensates for the reduction in c,,. 
_ The variation of c,,/2k,,, versus frequency for the various force levels, is 
ee. sfc shown in Fig. 11. This ratio increases with increasing force. Note that over es 
range of frequencies studied, the ape ‘nearly straight lines with 
a 
z 


Pim. 


. 


400 
My HORIZONTAL FORCE (hips) 


of _/2k, with Level of Force, for Pile Diameters a 
kip = 4.45 KN; 1ft= 0.305 m 


the lines i increases. 


a or offshore structures under seismic excitation, it would — 

that relatively simple formulae could be obtained from the results of this study. 

No attempt is made, however, to obtain such formulae, ‘since the main purpose 

of this: “study is to investigate trends which may not be affected too much > 

_ by the simplif ying assumptions, rather than to obtain specific quantitative results C2 


| | 

| 

| 
| | 


In the low frequency range, the radiation would be very small and 
the only | important contribution would come from the internal soil damping. 
This would be the range of interest for offshore structures under extreme wave ra 
_ The effect of the pile diameter in c,, is shown in Fig. 12 for three diameters — 7 
(d = 4 ft, 2 ft, and | ft), and for constant frequency 1 Hz. The variation 
of c,, is plotted versus level of force. Typical behavior, associated with small . 
PA frequencies of initially increasing and then decreasing c,,, with increasing level 
_ of force, is evident for the three diameters. The trends are very similar in | 
_ The variation of c mS 2k,,, with level of force for the three pile diameters, 
24 
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FIG. 14.—Variation « of k,.. and c,, /2k,, with Number ‘of Cycles, N (1 Ibt/ft = 


14.59 N/m; 1 kip = 4.45 in: ft = 0.305 m) 


is shown in Fig. 13. A monotonic increase with level of force | ise ; common to 


All the results considered previously were obtained assu assuming 10 cycles 
pe for the constitutive model described earlier. As was pointed out, the 
_ number of cycles will have some effect on the values obtained and the comparison 
with the p-y solution. To investigate the relative magnitude of this effect, the | 
“ft diam pile is studied under a force of amplitude, 600 kips (2, 670 kN), and 
“the ‘equivalent dashpot constant, and the effective ‘damping, V versus 
number of cycles, N,forthiscase, 
1 
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An attempt has been made in this study 1 to i investigate the effect of nonlinear ear 
Ed soil behavior on the dynamic stiffness of single piles, using a simplified nonlinear 
> model and an iterative linearization procedure. Because of the crudeness of be. 
7 the soil model, and the fact that important effects such as slippage, separation, s 
and gapping are not considered, care must be exercised in interpreting the results. *: ie 
oe should not, in particular, attribute too much significance to the actual 
"numerical values. It appears, however, that some important trends can te 
_ predicted well, even with this simple idealization, fa 
Comparing the results of this analysis with those obtained using the p-y curves — 
(perhaps the best solution available and with an experimental backup), it appears — . 
that the 1 main is due to the lack of tension in soil, 
in future studies to use a more . realistic soil model [such as the multiple- yield © 
‘hen idealization of Prevost (14)], and a three-dimensional solution, it is’ a 
% perhaps more important to reproduce in the model these other nonlinear effects. 
a The results indicate that as the level of forces i increases, there is an increase £ 
a in internal soil damping, but a ‘decrease in radiation. Still, when both effects 
are combined, the effective damping ratio tends to increase ( (while the — 
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Appenoix 


The following symbols are 
shear wave | velocity i in soil; 
imaginary part of horizontal pile stiffness dashpot co 
for horizontal excitation); ages of | Nolens 
eh . = imaginary part of vertical pile stiffness (equivalent dashpot constant 


/ 2k, = effective ratio at each frequency; 
soll | damping ratio; | 


| 
| 
hittin 


= = 
Fo = horizontal force at top of = 


= 


real part of horizontal pi e sti ness; 


real part of vertical pile stiffness; * dn 


p= soil resistance e (force per v ‘unit of pile); 


t = wall thickness of pile; goed bow q 


x distance from mudline; \o 
x, = depthofgapzone; 


Ye deflection at one-half > 

z = verticalaxiss MoT 

r* = shear strain amplitude; nd tho? Xe 


at one-half ma 


, = a stress of failure for static loadin - 


4 E, = 
at small strains; | 
— 
/ 
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CoupLep HEAT AND WATER F LOWS 


Omar N. Abdel-Ha Hadi,' A . ASCE, and James “ASCE 


The phenomenon of water - movement in soils due to thermal gradients 7 ; 


- important in relation to several problems, including heat dissipation from buried 

{ bry cables, underground storage of hot and cold materials, moisture movement 
and distribution in pavement bases and subgrades, and the ‘soil as a a medium 
for plant ‘growth. Thermally- “induced ‘moisture movements c can lead to ) changes: 

. ‘in both the thermal and isothermal properties of the soil which can subsequently — = 
affect the functioning of the soil for its intended purpose. __ 
_ Thermally induced soil water movement in unsaturated soils is a complex 
process. The various components of thermally-induced soil water movement > 
an the physical laws governing them are presented in this paper. The relationships 


developed are illustrated by their use as a basis for predictions of transient : fi 
nd steady-state temperature and moisture . distributions around buried electrical 
— 


transmission cables. aod DoW vies the followine sauauon 
- Buried transmission cables generate heat, causing their sheath temperature _ 


~ to rise. The allowable current loading may be limited by the maximum permissible “d 
cable sheath temperatures. T cmporateres higher than 50° C-60° C may lead o a 
cable failure and thermal instability | of the soil around the cable, caused by 4 
the fact that the soil water migrates away from the cable under thermal gradient, 
Jeading to a marked decrease in thermal conductivity of the cable environment. 

_ The result of this is that the ability of the backfill to dissipate heat is decreased. = 


The of. in unsaturated soils is less understood 


than moisture movement in saturated soil. This is due to the simpler a A 4 


"J mechanism and mathematical description associated with saturated flow as 
- opposed to unsaturated flow. In saturated soils moisture movement occurs i 
‘Geotechnical Engr., Geotechnical Consultants, Inc., San Francisco, Calif. 
-? Prof. and Chmn., Dept. of Civ. Engrs.» Univ. of California, Berkeley, Calif. 94720. a 
Note.—Discussion open until April 1, 1982. To extend the closing data one month, 
a written request must be filed with the Manager of Technical and Professional dese seg 
ASCE. Manuscript was submitted for review for possible publication on September 12, 
_ 1979. This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings 
of the American Society of Civil @ASCE, V Vol. 107, No. GTIl, November, 
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_ A comprehensive review of the literature on water migration under thermal hy 
gradients showed that several theories now exist. To predict the rate of water 
migration, many researchers have used Fick’s law modified to take into account 
_ both the reduction of diffusion cross section by the solid matrix and the liquid — 
_ water and the tortuosity of the diffusion path through the medium. Thus, va 


‘ the | liquid phase only. In unsaturated ‘soils t the flow ies place had we the 


a Taylor and Cary (18) used the general theory of icreversible wpe 
to account for the simultaneous flow of water and heat in porous media under 
7 the influence of thermal gradients and other secondary energy gradients. Philip Bf 
and DeVries (13) considered moisture flow in vapor and liquid phases, and 
the interaction of liquid and vapor phases. Laboratory experiments, Casset, 
5 et al. Q), showed close correlation with the Philip and DeVries theory, ‘while 
, Taylor and Cary’ s approach underestimated the flow rate by a factor of 10 
to 40. The theory of Philip and DeVries was modified and used for the analysis — 
_ of moisture flow in base courses under pavements by Raudkivi and U’u (14). 
- 2 theory of Philip and DeVries appears to provide the most Le 
basis for development of a model for prediction | of _ transient flow of 
in unsaturated soils under thermal gradients. “Accordingly, ‘it ‘was adapted for 


the tudies described h 


In this paper, a summary of the Philip and DeVries theory is presented. 
me of the soil parameters needed, according to the theory, is described | 
along with examples. The relative importance of the properties affecting the 
rate and amount of water migration under thermal gradients is illustrated. Finally, ; 
a field test program is described which | provides | a basis for comparison between = 
-_ predicted and actual values of moisture and temperature distributions. = 
In an analysis of moisture movement in unsaturated soils under a thermal — 
- gradient, the following must be considered: (1) Vapor phase transport; (2) liquid 
phase t transport; and (3) interrelationship of the vapor and liquid phases. pris ae) 


G v apor Phase Transport. —The following equi equation was s obtained by Philip _— 


( _ modified so as to ‘apply to porous media, and extended to separate the ane 
and thermal components of vapor transfer and to take into account the effect 


in which: Qvap = Vapor flux ails in grams per square millimeter per second; — 


— * = density of water, in grams: per cubic millimeters; and Dr, = thermal | 


however, that this theory underestimates the quantity of water vapor transport _ 
| 
= 
| 
4 
| | 
| 


ao vapor diffusivity, in square millimeters per second per degrees Celsius, which — 


4 g in which D, = molecular diffusivity of water vapor in air, in square millimeters q 


per second; v = “mass flow factor’ introduced to allow for the ‘mass flow — 
of vapor arising from the difference in boundary conditions governing the ai g 
of vap g y g the air 
and vapor components of the diffusion system. It is equal to P/(P — p); in 
which P = the total gas pressure; and p = the. partial vapor pressure in the © 
- pore space; a= tortuosity factor allowing for extra path length; a = volumetric 
- air content of the + soil, in cubic millimeters of air per cubic millimeters; h = 
‘relative humidity; Po = density of saturated water vapor, in grams per cubic 
millimeters; = temperature, in degrees Kelvin; V = partial derivative. with 
‘respect to distance; and D, = isothermal vapor diffusivity, in i in square millimeters. 
Pet second which is equal to to = ati 


T 96 


R= gas const, in ergs per grams per Celsius; = = = suction head, 
in millimeters, in thermodynamic equilibrium with the water in the medium; 
= _ @ = volumetric water content, in cubic millimeters of water per cubic millimeter 
- of soil volume = w Pa/ Pui w = weight of water / weight of dry soil, as a percentage; = 
p, = soil dry density, in grams percubic millimeter, $= 
iF Liquid Phase Transport. —Philip and DeVries obtained the following equation - 


_ for the liquid flux density. This equation is based < on Darcy’s law for liquid 


_ transfer in unsaturated media ant —— to separate e the isothermal and thermal — — 


in which iia iq = liquid flux density, in grams per square millimeters per Oe 7 
- D, = thermal liquid diffusivity, in square millimeters per second per degrees 


ke unsaturated hy ydraulic conductivity, in millimeters per second; 
o = surface tension of water, in dynes per millimeter; D, = isothermal liquid s 


in square second, which is also equal to 
‘in which pis in the positive z direction. « “> 
___ Interaction of Vapor and Liquid F Phases. - —Philip and DeVries (13) modifi ‘= 


Ulu 
| 
| 
4 
ne equation for the thermal vapor ivity T, given ne Simple Vapo : | F 


bys a factor of 3 to 10, as mentioned earlier in this pa paper. per. Their modified equa equatio: 
accounts for the interaction between vapor and liquid phases, and the difference 
between average temperature gradient in the air-filled pores and in the soil 
as a whole. They found the modified equation to predict order of magnitude 
general behavior in satisfactory with experimental o observations. 
me it 


Th fi 


square millimeters per second per degrees Celsius. ...... 


which f(a) = a/a, for 0 < = for a = a,; in ‘which a, = the 


value of a (volumetric at which liquid continuity in pores | no longer a 
exists; or a, = the value of a at which 6 = 6,; in which 6, = the moisture > 
"content at which k (hydraulic conductivity) falls to some small ‘<stitreny fraction — 
Bed its saturated value, and b = the ratio of the average temperature gradient 
in the air-filled pores to the overall temperature gradient. Ra 
Comprehensive Theory.—The total water movement in an unsaturated = 


_ due to a temperature gradient and its resulting water content gradient is equal 
: to the sum of the flows occurring in the two phases; vapor and liquid. Thus: _ 


in which g q= the total water flow, in grams per square millimeter pe per second. 


+ Dy, = = water diffusivity . (11) 


Eq.9 9 is the oration expression for moisture movement under a thermal gradient 
in unsaturated soils as proposed by Philip ‘and DeVries. Differentiation of (9) 
and application of the continuity requirement gives the 

a 


The heat conduction equation for the soil is: 
in which cs = volumetric heat capacity the soil, in per lite millimeter 
conductivity, in calories per millimeter p per 


= 
| 
j 
in which Pw «2D Theemel eater 0) 
D, =D. 


per 12: and 13 a are the ‘describing coupled 
heat and moisture transfer in partially saturated soils. 8” or" a 
Evatuation or Son Parameters od} Yo ewibar 
- From Eqs. 12 and 13 it is clear that the variation of four soil parameters 4 
with the water content (or with degree of saturation) is needed in order to 
predict the moisture transfer. These parameters are: (1) k, (hydraulic conducti- — 
vity); (2) ; D, (isothermal water diffusivity); GB) Dr (thermal water diffusivity); 4 
and (4) x - (thermal diffusivity). These four parameters have been evaluated for 7 
several soils (9) as well as for the soils of the field test described later. Crushed 
limestone, referred to as PG and E limestone screenings, were used in the 
field test, and the soil parameters obtained for this soil are eel baal 
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ee range used in the evaluation of the soil parameters was 10° C-60° c. 
The methods for evaluation of the soil parameters | are given in the following 04 
Hydraulic Conductivity, k,.—Marshall (8) derived the following equation for 
‘the relation between the hydraulic conductivity and the size distribution of the 


in n which k = specific hydraulic conductivity, in square millimeters; « = porosity, — 


in cubic millimeters per cubic millimeters; n = total number of pore classes; > 
and r, = mean 1 radius of the rores in the pore class i, in millimeters. Eq. 


| 
_4 
x 


> 
14 it possible to the hydraulic of potous 
from the size distribution of its pores. Pore size measurements can be made 
_ by measuring the water withdrawn when the suction on the soil is ——- 7 


increased. radius of the largest water-filled a suction of 


whi o = surface tension of water, in per Pw = density 


on Head, (cm of woter) 


cho 


> 


Su 


ms 


vag 


FIG. 2.—Relationship Between Suction Head and Volumetric Water Ci Content for PG 
&E Crushed Limestone 


of water, in grams per cubic millimeter; g= = ot due to gay, in 


a millimeters per square second. In applying Eq. 14, it is often more convenient 
14 becomes 


7 _ to use the suction in place of the pore radius. Therefore, Eq. 
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The units of the hydraulic can converted to millimeters per 


second by multiplying the hydraulic conductivity by 7 in which 
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derivation of this by} Marshall, as the filled 
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FIG. 3.—Variation of Hydraulic Conductivity of PG & E Crushed Limestone with | 
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porosity at water content 0, for which hydraulic conductivity (kb is 
calculated, and n was the number of pore classes in the water content interval _ 
zero to 8,. Thus both and n decrease as the water content ied 4 

. Volumetric Woter Content, 8 oni 
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Unit OF KN/m (US pct) — 
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7 - Green and Corey (5) proposed another approach which is similar to Marshall’s 
and leads to the same results. They assumed € to be the water-saturated porosity, 


‘€ and n to be the total eunter « pore classes at « = 6,. In this ‘uppreach 
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lic is for 


“different values of water | ‘content (8). This approach was used in calculating 
tthe: hydraulic conductivity of the soils in this investigation. 
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FIG. 5.—Variation of Isothermal Vapor _Diffusivity of P of PG & & E Crushed Lir Limestone a 

_ A matching factor is usually introduced into Eq. 17 to equate the calculated 


3 and measured hydraulic conductivities. Elzeftawy and Marshall (4) stated that : 


4 t wa. 
| 
— 
we 


NOVEMBER 
; matching s at water saturation has a adistinct advantage over matching at ¢ at desaturation _ 


water contents mainly because it is simpler and more accurate to measure saturated _ 
hydraulic conductivity. Rewriting Eq. 17 in the a form and introducing 
2b, is =1, 
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FIG. 6. —Variation of Thermal Water Dittusivity of P PG & Crushed with 
Degree of Saturation pane 


in millimeters p per second; i = the last wk content class on ~~ wet end, 
" e.g., i = 1 denotes the pore class corresponding, to the saturated water content, os 
and i = m denotes the pore class corresponding to the lowest water content 
for which hydraulic conductivity is calculated; (k,/k,.) = matching factor = 
_ (measured saturated hydraulic conductivity /calculated saturated ma con- a 


u ductivity; and m = total number of pore classes (pore class is a pore size 

Tange ‘corresponding to a water | content increment) between 6 = 6, and 6,, 

in which 6, = the lowest water content on the experimental moisture content-suc- 
tion head curve Th hus: ¥ 
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6s wit at ai ayats 


Constant value for m is used at all water contents, and the value of m establishes a 


types). wy value of m = 12 was en Rare to calculate the results given in this 
_ From the comparison between experimental and calculated values, it has been _ 
GB, 5, 7, 11) the matching» factor, 
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FAG. 7. —Variation of Thermal Resistivity of PG & E Crushed Limestone with Degree 
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successfully predicts the hydraulic conductivity for a wide range of 7 a me 
The suction head-water content relationship for sandy soils can be obtained - 
wang ‘Tempe”” "pressure cells. In this device cylindrical s: samples approx 50 


- ‘atmosphere. After the soil is placed in the cell, it is saturated with water. 
The soil is allowed to drain following sequential subjection to increasing air 
pressure. Volumetric water content is determined from the value of the withdrawn © 
- water after equilibrium is reached at each pressure and the water content of 

F the sample at the end of the test. The suction head-water content relationship — 

- for PG and E crushed limestone is shown in Fig. 2. The suction head- water 

content relationship is a hysteretic function; i.e., 8 depends upon whether the 


gg at saturation. m should be chosen > so that the water content-suction head 
+ 
| 
= 
of Coble Temoerature on Volume vin Werer 
¢ 
| 


soil is drying or wetting. The data in Fig. 2is ~~ drying. ee 
iD ‘The suction head-water content curves are strongly affected by the soil gradation | 
and the soil fines content (silt and clay). For silty and clayey soils, the mene 
_ range of interest extends to much higher values than can be measured by the 
“‘tempe”’ cell. Therefore, for finer soils suction-water content measurements 
_ were made using a pressure plate apparatus which can accommodate suctions 
as as 15 atmospheres = a 


_ was measured for the sandy soils in the laboratory 1 using falling head permeability 
test equipment. Eq. 18 was used to calculate the hydraulic conductivity using — 
a short 7 program written for this purpose. The hydraulic conductivity 
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FIG. 8. of Thermal Dittusivity of PG&E Us 
of PG and E crushed limestone is shown in Fig. 3 as a function of both water = 
Isothermal Liquid Diffusivity, D, o-— Eq. 6 was used to calculate D, —_ 
the soil water content suction head relationship is hysteretic, it follows that 
‘D, is also an hysteretic function. Then, in calculating D, from k, data (Fig. 
3), it is necessary to consider whether the soil is wetting or drying. The values 
of D,,, for drying, versus + pee) of saturation of PG and E E crushed apace 
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Isothermal Diffusivity, D, was calculated 3. 
The D,, versus degree of saturation relationship for PG and E crushed limestone fee 
is shown i in Fig. 5. The value of D, ,, Must decrease very sharply at low degrees — 
of saturation, because the relative humidity (A in Eq. 3) from more 
n zero between S = 10% and S = 
Thermal Water Diffusivity, D, 
Thermal Liquid Diffusivity, D;,.—The D, was using Eq 5. The 


surface tension coefficient {1 /a(do/dT)| can be taken as a constant 
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Vapor Diffusivity, D,,.—The D,, was using Eq. 7. 
thermal water diffusivity, D,, ‘versus wy hin of saturation for PG and E crushed 4 
limestone is shown in Fig. 6. ‘The both the thermal liquid diffusivity, 


Eq. 5 we can see ‘that D, vert. ye cena zero at S : = -0 and 100%, because 
_ the hydraulic conductivity, k,, is practically zero at S = 0, and the suction — 
head is zero at S = 100%. From Eq. 7 it is clear that D,, approaches zero 
= 0 and 100% because the relative humidity, A, "decreases to zero at 
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ian” yi S = 0, and the volumetric air content (a) and the function fa) a are eal to 


x=— 


A, was evaluated ‘in the lat laboratory using g the pam 4 
needle method (10). The thermal resistivity (1/A) ore ‘Cages of saturation | 
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FIG. 12. Variation of Volumetric W Water Content Pro} Profile with Time with 


for PG and E crushed limestone is shown i Fig. 


_ The vi rr heat capacity, C, can be determined from the expression: a 


calories per cubic millimeter per d Celsius 


in which c= specific heat of dry ool, in calories per gram; y, = unit Pom 
= water conte ercent dry 
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thermal versus degree of saturation relationship for PG and 


_E crushed limestone is shown in Fig. 8. The shape of the thermal diffusivity — 
- (d-degree of saturation (S) curve reflects the changes in the ratio A/C with | 
increasing Ss: saturation. The values of \ increase rapidly at low degrees of saturation © 
(less than about 30%), while the rate of i increase of ‘decreases appreciably 
2 at higher saturations. On the other hand, C increases steadily with the increase | 
_ in the degree of saturation. Therefore, the values of \/C increase at low mee 
a na y of saturation, until a limit is reached at which the ratio A/C starts to decrease. | 
pie _ The preceding relationships and properties make possible the prediction of . 
Moisture movements under combined | hydraulic and thermal | gradients. Their 
is. liable has been investigated as a part of a research program on the 
Be dissipation of heat from a buried electrical transmission cables (9). The complexity a 
of the relationships and requires the use of puter 


‘qe 


“AIG. 15 —Deson of Moir and and Temperature Sensors in Field T Trench 
a Gibines. An explicit finite difference computer program was developed for | 
* nalysis of flows in one dimension (1). Although a two-dimensional program — 


rig i. ould better r represent the e actual conditions, its greater - complexity prohibited | 


its at this stage of the research. 


“The ‘aula obtained for three soils which have been widely used as cable ia 
backfill materials (PG and E crushed limestone, Round- Robin sand, and Fire 
Valley thermal sand) were used to analyze the relative importance of different — 

- parameters on the amount and rate of water movement under thermal gradients 
ach x trench, cable, , and —— system using =e one- ee coupled flow 
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_ which the cable and backfill are placed, were > taken as 4 ft of highly hed 
silty clay overlying a l- -ft thick layer | of low plasticity silty clay underlain by 
silty sand and gravel. This profile is essentially that at the field test site to 
be described later, and it is referred to as the RFS soil profile. == = 


FIG. 16. —Measured and Predicted Temperature Distribution 5 Days After Heating ; 
The following parameters were studied: 
1. Cable temperature. 
2. Backfill soil type. 
3. Open or closed system at ground surface.  < 
Depth of watertable. 20 bess 
6. Values and distribution of initial soil water content. 
. Values and of initial soil temperature. ath 
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--9. Weather variables: (1) Air temperature; (2) wind ae and (3) ‘solar 
a Only the effects of cable temperature and backfill soil type are presented — 
here to illustrate the type of results obtained. The ground surface was assumed _ 
to be at constant temperature and open to free evaporation. It was assumed © 

_ also that the cable sheath temperature was increased instantaneously toaconstant _ 
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| «FIG. 17. —Measured and Predicted Temperature Distribution 20 Days After Heating 

9 | Shows the predicted change of v volumetric water content distribution 
with time for different cable temperatures and for the conditions shown #7 
the sketch drawn beside the results, and Fig. 10 shows the change in moisture © 
- content of the soil adjacent to the cable with time for the same conditions. — :/ 
From Fig. 9 it is clear that the moisture content of both RFS soil layers do 
not significantly either with time or with is also | clear 


| 
| 


a os ones its initial as-compacted value with time for all cable tempera- a 
y tures. This is mainly due to gravitational flow, because the soil cannot hold fl 
_ more moisture than its field capacity (maximum amount of moisture that can _ 
be held against the force of gravity). From Fig. 10 it can be seen that after — 

100 days of operation, the moisture contents of soil adjacent to the 60°C and © 

90° C cables decrease to ; . about 62% and 40% of that around the 15°C cable, 

respectively. Such moisture losses could be important if they are accompanied — 

by decreases in thermal conductivity. In that case the ability of a cable backfill — 

system to dissipate heat could be impaired. ae 
Figs. 11 and 12 show the effect of two different backfill soil types [PG 


; , and E soil and RFS soil (2 ft)] on change of moisture content profile with 


time for the. conditions shown. The values of the parameters needed for the 
predictions for the RFS soil are presented by Abdel-Hadi (1). Although RFS 
soil can hold more moisture than PG and E soil, it would not be a good backfill 
_ material around buried cables for two reasons: (1) It is a silty clay having oe 
relatively high thermal and (2) high densities are difficult to obtain 


The results show that the cable temperature and soil type exert the greatest 


influence on moisture migration. The cable depth may have a relatively large — 


effect for an open system at ground surface (6 = 0 at ¢ > 0). The type a 
—- at ground surface (open or closed), water table depth, initial ground | 
. _ temperature, and initial values and distribution of the volumetric water content _ 
4 “do not have significant effect for the conditions considered for these analyses. 


In general, short term weather variations should not have important effects. 


: on the moisture changes around the cable. — 


A 


a theoretical and actual heat and w water flows ‘around buried power ‘cables. The 


a test site is located at the University of California Richmond Field Station (RFS). 
The trench analyzed is one among ten trenches constructed for analysis of — 


r “the effectiveness of different backfill treatments to enhance heat dissipation | 


and to evaluate the accuracy of the developed heat and \ water flow predictive _ 


~ Five bore holes were drilled in the test site to determine the soil conditions, . 


and laboratory tests were performed on the soil samples obtained to determine _ 

4 the in-situ moisture contents, densities, and soil classification. A typical boring 
jog is shown in Fig. 13. The soil at the test site consists generally of a stiff — 

dark brown, highly plastic silty clay (CH) overlying a slightly inclined layer 
of light brown low plasticity silty clay (CL) about | ft thick at a depth of a 
about 4 ft. Below this layer is a layer of silty sand and silty gravel with some 7 7 g 
- silty clay seams. Free ground water was encountered at a depth of about 12 

_ Thermal and hydraulic properties of the three main layers of the test site 

— evaluated in exactly the same way as described earlier for the PG and 

E crushed limestone and are reported by Abdel-Hadi (1). mort 

: Conduit and Heating System.—An 11-ft long, 8-in. 1.D., carbon steel standard 
7 pipe was used as the casing of a simulated electrical transmission cable. A 
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heating: rod was, asa ‘element inside ‘the pipe to the heat 
that could be expected in a real buried cable system. . The space between the 
heating element and the pipe was filled with oil. Fig. 14 shows a schematic 
—- of this heating system. The system is capable of producing 200 watts /ft 
under full load. A cont controller was designed to maintain constant temperature 
on the pipe, supply constant heat in input, or supply | programmed heat input. odT 
a *y avoid the build-up of high pressure due to thermal expansion of the oil, 
én oil tank with a relief valve (+10 psi) was installed near the ground surface. 
_ Weather Parameter Measurements.—The following weather parameters were 


Wind speed and direction. a an the 
3. Net solar radiation. 


Precipitation. 

In-situ Measurements of Soil and Backfill Water Content.— Accurate moasuse- 
‘ment of water content of soils is essential for studying the problem of water 


q movement in- ~situ. A Teview of the methods tl that could be used for this purpose 


is based on n the fact that the dielectric constant, which determines the capacitance, 

depends on the water content of the soil. The dielectric constant of water is 

about 80 and the dielectric constant of oven-dry soils is only about 2. 6. Thus, 
the. addition of water to soils causes a significant increase in | the ‘dielectric: 


. i g 16) and Wobschell ( (20) were used. A typical probe. consists of a Sasblia 7 
case 25.4 mm (1 in.) diam by 152.4 mm (6 in.) long containing an electronic 
circuit inside. The probes are calibrated in the laboratory using the same soil 

be encountered in the field 
‘Field Performance of Moisture Probes. —Some of the ‘moisture Probes failed — 
ai due to leakage of water into the electronic circuit inside the probe; i.e., ‘insufficient 

4 water proofing, and some other probes failed or gave misleading results due ! 

- to temperature effects. Modified moisture probes have been used in subsequent 

trenches, with better water proofing and with components in the electronic 
circuit rated to operate satisfactorily at temperatures as high as 120° Sapper 
Construction of Field Test Trench.—The field test trench was constructed o 


' on May 2, 1978. It is 0.61 m (2 ft) wide by I. 22 m (4 ft) deep by 3.05 m i 


(10 ft) long and w was backfilled with PG and E crushed limestone s screenings. 

4 The screenings were placed at 18.1 kN/m?* —18.8 kN/m? (115 pef-120 pef) dry | 
density and about 13% water content. The instrumentation around the cable . 
-: > consisted of 25 temperature sensors ae 10 moisture probes arranged as shown © 
in Fig. 15. Measurements of the soil and weather data were initiated immediately : 
after the construction of the trench. Measurements under isothermal 


continued for 20 days. Then the pipe was heated to a sheath h temperature of 


60° C, and this temperature was held constant. 
Comparison of Measurep aNd PREDICTED TEMPERATURE 0} sub 


Eye. 29 shows tbe values of theoretical 
The purposes of the field test program included of theoretical 


methods for or predicting temperature dist distributions as af of time and Position 
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““HEAT”’ (17). This program solves both transient and steady-state linear heat j 
conduction problems in two dimensions with time dependent boundary conditions. a 

The governing equations in the program.are the heat balance and heat flow. — 
The effects of conduction, convection, and solar radiation at the g ground surface 
% can be accounted for. Predictions were made for the time required for the 
: pipe sheath to reach 60° C and for the temperature distribution around the pipe 


4 


FIG. ‘18.—Measured and Predicted Distribution 40 Days After 
at different times after | pipe heating was ‘initiated. Measured weather data a 
P - initial moisture and temperature conditions in the field were used in the predictions. ; 
The predicted time required for the cable sheath temperature to reach 60° . 
due to a heat input rate of 200° watts /ft was” about 11 hr, while the actual 2) 
_ Figs. 16-19 show comparisons between measured and predicted temperature 
distributions at 5 days, 20 days, 40 days, and 80 ead after heating | the pipe, 


/ 
/ 


Also ‘shown on each figure a are the 
temperature, V,, = wind velocity, J = solar radiation) and soil thermal resistivities, 
p, used for the awide Sart formed by NET WO 

_ From these figures the following can be ol observed: rt 


or times less than 40 days after heating the cable Sea Jalatniais is me 
between the measured and predicted temperature contours. 

ee For times greater than 40 days after heating the c: cable, the ag agreement is 

, but the | gap between the measured and predicted contours 
increases with time. This is mainly because of the accumulating error resulting 
from using the temperatures obtained theoretically at any time as the initial 
value for prediction of the temperatures at the following time. “eat aa 


t In general, the predictions of the temperature distributions around the cable 
= program HEAT are good. The temperature profiles generated by the 
one-dimensional computer program for the conditions presented in Figs. 16- 19 
a _ match fairly v well ll with the 1 measured field temperature profiles above the cable. 
the | program underestimates the temperature gradients measured 
the trench immediately below the cable. This is believed to be de due to tl to the 
two-dimensional nature of the problem. | 


Comparison of Measureo ano Previcteo Moisture Content Distrieunions 


the soil expected to be used as arty After the heating of the field trench, oT 
the results obtained from the moisture probes indi agers that t the drying of ‘the 


_ than predicted. Several possible reasons for this difference i in ‘iocea rate were uf 

_ The saturated hydraulic conductivity of a ouuate obtained from the actual 

_ backfill soil, Sample 2, was found to be | x 10~* mm/s compared to 2. a4 


: x 10~* mm/s for the sample of the soil tested previously in the laboratory, — 


: Sample 3 Also, Sample 2 had slightly “more fines p passing #200 sieve (about 
> 8%) than Sample | tested previously (about 5%), although the “two ss samples 
_ had essentially the same gradation curve. These differences are very important, 
because the predicted values of water content are very sensitive to the hydraulic — 
conductivity of the soil. Using the new value and a new measured relationship 
between the suction head and water ¢ content for Sample 2, the coupled flow . 
parameters were recalculated and used to - predict new values for the water 
, = Because of the failure of some of the moisture probes due to both insufficient | 
-_- waterproofing and temperature effects, and because of some suspicious water 
' content values obtained from some other moisture probes, a borehole was made 


every 20 days to obtain actual water contents from soil samples. Bm 


Fig. 20 shows the measured values obtained from the bore hole oie 
a and the predicted values using properties weagened from Samples | and 
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FIG. 19. —Measured and Predicted Temperature Distribution 80 Days After Heating © 
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- 20. wey of wel roth and Predicted Moisture Content of Backfill Soil 


os measured profile lies inside a a wide band fi ormed by the two predicted profiles. " 
2. Neither of the predicted water content profiles shows the drop in water — 


content just above the cable (water migration under thermal gradient) which - 


 — ae relatively small change in the soil gradation and density may have great 
effect on the predicted water contents. This fact raises questions about | 
ability to predict accurately the amount of water content changes under thermal os 
gradients for materials of the type studied herein, especially considering the 
uncertainties involved in both measuring the properties of truly representative — 
samples and the field handling = placement of soils. Further research on — 
this « question is in Progress. Klemen for 
Dimensionat Versus Two Dimensionat Move. 3 


he geometry of the trench-cable system allows the flow of water and heat 


pers from the cable in all radial directions. Nonetheless, it had been thought, 


since the cable ‘is confined in a trench surrounded by relatively impervious 


soil, that some insight into the moisture. migration | problem and some estimate 
a. water content profiles in the test trenches might be obtainable from a 
one-dimensional model. Additionally, a one-dimensional formulation was used 
first because of the greater difficulty associated with the two-dimensional — 
‘numerical analysis. Final c conclusions about the e adequacy o of the ¢ one-dimensional — ; 
- ‘model cannot yet be « drawn from the analysis of this single tr trench. ole il 4 
Development of a two-dimensional model would provide the following advan-_ 


tages as compared with the presently available one-dimensional solution: = =—> 


. Actual geometry could be simulated. bee 


2. More representative boundary conditions could be used. len A 


3. Comparison between measured and predicted water content and temperature 


‘@ The theory of of Philip an and anaes 013) was used for prediction of of thermally 
; a induced moisture movement in unsaturated soils. Theoretical analysis s of moisture 
movement and the parameters affecting it in a typical trench, buried electrical e 
ie cable, and backfill system indicated that the most important parameters affecting _ a] 
the rate and amount of moisture migration away from the cable orang be 
cable temperature and soil type. Gravitational flow, surface evaporation, and — 
a, isothermal diffusion have great effect on the moisture changes i in the backfill 7 


field test an actual buried cable system was conducted to evaluate 


_Capacitance- probes for moisture content measurements | did not perform well “i 
beens of insufficient waterproofing ; and temperature re effects. Actual moisture — 


Neither of the predicted water content profiles, using properties of Samples 
| | 
| 
&g 

and are PeSDOUSIDIC TOT Caliv Chanees Of ihe Moisture if ine 


NOV 
from two backfill differing slightly in gradation and 
_ matched poorly with the measured profile. However, the measured profile fell = 
_ within the band formed by the two predicted profiles. ‘The predicted water 
values were found to be very sensitive to small “changes | in the soil 


The finite element computer program ‘“‘HEAT’’ was used to predict the 
temperature distribution around the cable at different times after initiation = 
cable heating. The results showed that there is excellent agreement between 
om measured and predicted temperature contours up to 40 days and good 
agreement between the contours from 40 days- 100 days. 

The studies described in this paper were sponsored by the Electric — 
Research Institute, Palo Alto, California, under the agreement No. RP 7841-1, 
Backfill Materials for Underground Power Cables. John C. McMillan a wadete 
‘Tesearch. of he ¢ University of California, at Berkeley, California, reviewed 


Abdel- Hadi, O. of Heat and Water Around Underground Power Cables,” 
_ thesis presented to the University of California, at Berkeley, Calif., in 1978, in partial . 
7 fulfillment of the requirements for the degree of Doctor of Philosophy. nara gers 7 
2. Cassel, D. K., Nielsen, D. R., and Biggar, J. W., “Soil Water Movement in n Response — 
& _ to Imposed Temperature Gradients,”’ Soil Science Society of America, Proc. 33, 1969, x 
3 Elzeftawy, A., and Dempsey, “Prediction Model for Hydraulic Conductivity 
and Moisture Diffusivity of Highway Soils,”’ presented at the January, 1977, 56th 
Annual Meeting of the Transportation Research Board. = 
4. A., and Mansell, R. S., ‘Hydraulic Conductivity Calculations for Unsatu- ih 
rated Steady-State and Transient-State Flow in Sand,”’ Soil Science Society merica, 
; 5. Green, R. E., and Corey, J. C., “Calculation of Hydraulic Conductivity: A Further 
Evaluation of Some Predictive — Science Society of America, Proc. 
. Gurr, C. G., Marshall, T. J., and Hutton, “Movement of Water in Soil Due 
to a Temperature Gradient,’ Soil Science, Vol. 74, 1952, pp. 335-345. _ Rae 
2 7. Kunze, R. J., Vehera, G., and Graham, K., ‘‘Factors Important in the Calculation — 
_ Of Hydraulic Conductivity,”’ Soil Science Society of America, Proc. 32, 1968, p -— 
a 8. Marshall, T. A **A Relation Between Permeability and Size Distribution of Pores,” 


Power Cables,’’ Phase I, Thermal Resistivity Measurement Methods, Backfill Treat- 
é ments, Heat and Moisture Flow Analysis, Report prepared for the Electric Power 
_ 10. Mitchell, J. K., and Kao, T-C. , “Measurement of Soil Thermal Resistivity,’’ Journal 
— of the Geotechnical Engineering Division, _ ASCE, Vol. (104, No. GT10, Proc. Paper 
ll. Nielsen, D. R., Kirkham, D., and Perrier, E. R., “Soil | Capillary Conductivity: 

_ Comparison of Measured and Calculated Values,’’ Soil Science Society of America, — 


Proc. (24, 1960, pp. 157-160. W fap 


rovided valuable tec 
| 
| 


Penman, H.L., 
: Science, Vol. 30, 1940, pp. 437-462. 
3. Philip, J. R., and DeVries, D. A., “Moisture Movement in Porous Materials Under — 
‘Temperature Gradients,” Transactions, American Geophysical Union, Vol. 38, The __ 
American Geographical Union of the Academy of Sciences, National Research Council, 7 
. Raudkivi, A. J., and U’u, N. V., ‘Soil Moisture Movement by Temperature Gradient,” 
Journal of the Geotechnical Engineering Division, ASCE, Vol. 102, No. GT12, Proc. 
Paper 12645, Dec., 1976, PP. 1225-1244. 


15. Rollins, R. L., Spangler, M. G., and Kirkham, D., “Movement of ‘Soil Moisture 
a Thermal Highway Research Board Proceedings, 33, 1954, pp. 
16. Selig, E. T., ‘Wobschall, D. C., Mansukhani, S., and Motiwala, A., Sy oe 
Sensor for Soil Moisture Measurement,”’ "" Transportation Research Record 532, 1975, 
ses ase critical part of (be. aciectian and design cf 
Taylor, R. L., “‘HEAT, A Finite Element Computer | Program for Heat- 
-_Analysis,”” Report 75-1, prepared for Civil Eng. Lab., Naval Construction Battalion — 
Port Hueneme, California, Department of Civil Engineering, University of 
18. Taylor, S. A., and Cary, J. W., ‘Linear Equations for the Simultaneous Flow of 
on Matter and Energy in a Continuous Soil System,”’ Soil Science Society of America, 


19. Taylor, S. A., and Cavazza, ae “The Movement of Soil Moisture in Response 10 _ 
_ Temperature Gradients,’’ Soil Science Society of America, Proc. 18, 1954, pp. 351-358. 
20. Wobschall, D., “‘A Frequency Shift Dielectric Soil Moisture Sensor,” JEEE — 


tions on Geoscience Vol. 16, Apr., However, 
of concern for many types ¢ nd 
determination of the setilemeat that 
members or architecture! elements, or both, is 
is affected by many factor, ihe type and sing 
an the properties of af the strectural me erialy ane the 
Ghiveria Tot tolerable have estabished the 
of of settionwn( and damage m esting baiidings; Usually, 


. be recugniaed. Te. primary problem the scarcity, af compleie 


pelichle cxse studies The sicomment of the type and degree of 
i deveiied of the and 


| 
| 


gut, pe 


3 


_ ame 


& 
| 
2 


JILDINGS 


of the enpertnces with these 


Settlement : analy ses are a critical part of the selection and design of building a al 
_ foundations. The anticipated settlement must be predicted, and, then, the tolerance 
of the structure to the predicted settlement must be evaluated. The purpose 
of this paper is to review current concepts and practices fi for establishing tolerable — 
settlements for buildings. Procedures for “prediction of settlement are not 


_ Previous reviews (3,4) have indicated many y of the difficulties of establishing - 
allowable settlements for buildings. Clearly, settlement that affects the safety 
or function of a structure is unacceptable. Such considerations may be the 
: governing factors for warehouses and industrial mill buildings. “However, — 
- appearance also is of concern for many types of buildings, and s significant cracking — 
a of architectural elements usually will be unacceptable to most owners: of 
7 _ The determination of the settlement that will cause significant cracking of 
_ ‘structural members < or architectural elements, or both, is a complex lesemiante 
analytical problem. It is affected by many factors, including the type and size — 
(of the structure, the properties of the structural materials and the subsurface 
7 ‘soils, , and the rate and uniformity of settlement. Because of these ‘complexities, — 7 
_ ¢fitical settlements have not been determined analytically. Instead, almost all 


- criteria for tolerable settlement have been established empirically on the basis a 


qualitative assessment of the level o of damage is correlated with some measure 
of differential settlement. Then, an attempt | made to relate” the 


@ observations of settlement and damage in existing buildings. Usually, a 
‘settlement to the maximum absolute settlement. 


a “Presented at the October 27-31, 1980, ASCE Annual Convention and id Exposition, held 


at Hollywood, Fla. 
"Prof. of Civ. Engrg., North Carolina State Univ., Raleigh, 


Note.—Discussion open until April 1, 1982. To extend the closing date one month, 
is written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on March 12, 1981. 
This paper is part of the Journal of the Geotechnical Engineering Division, Proceedings — 7 
of the American Society of Civil Engineers, ©ASCE, Vol. 107, No. GTll, November, — 
ISSN 0093-6405/81/0011-1489/$01.00, 
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everal inherent limitations to this empirical approach to tolerable settlements 
—— puld be recognized. The primary problem is the scarcity of complete and 7 
; g able case studies. The assessment of the type and degree of damage usually ; oA : 
— is qualitative. Often, detailed descriptions of the structures and soil conditions Fi . 


are not available. Sometimes erroneous assumptions are made regarding the 
relation of damage to settlement. Not all cracking is due to settlement, and, 
conversely, the absence of cracking does ‘not necessarily mean there has been 
no settlement. Relatively few settlement records are reported for undamaged 
buildings, and, consequently, the available case studies are or toward 
situations in which damage has been observed. 
With the preceding limitations in mind, the currently accepted settlement criteria _ 
il be reviewed and the basis for these criteria will be discussed. Some simple _ 
: conceptual ideas will be used to illustrate the | 


to aid the ae of observed data. 


effects of various factors re and 
Derinmions agiest bar nol reise of) to © oun 


ettlement to the vertical displacement of the 


lp 


moo 1. —Definitions of Settlement Terminology bivode 
i ei commonly is related to variations in the settlement at different locations 


in the building, the differential settlement is usually of more interest than absolute 
settlement. Common terminology for describing absolute and differential settle- 
ment is defined for cases without and with a component - of rigid body rotation - 
7 or tilting in Fig. l(a) and 1(b), respectively. The following terms are defined _ 
= vertical displacement at i; 8, = differential settlement between i and j; _ 
a A = relative deflection = maximum displacement from a \ straight line connecting be 
two tilt = rotation; By = = y- w= = angular 
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“ practices in Chicago. The period between 1880 and 1900 is of particular interest — x 
B.2 because of the common use of isolated spread footings and the recognition 
Of settlement as a major design consideration. Peck identifies 60 buildings of 

fe nine or more stories, which were constructed on spread footings during the | 


period | 1883-1895. Three major structures (the Auditorium, the Board of Trade, 
and the Polk Street Station) quickly ‘suffered major damage because of large 
differential settlements between heavily loaded tower sections and adjacent lower 
_ sections. As a result of the experiences with these structures, engineers rend 
"quently attempted to design structures with loads more evenly distributed so 
4 as to produce relatively uniform settlement. Then, the foundations were con- 
structed above the desired final grade to provide an allowance for settlement. 
For several buildings, including the well-known 16- -story - Monadnock Building 
and the 21l-story Masonic Temple, allowance was made for 8-9 in. (200-300 # 
of anticipated settlement. 
Peck summarizes that as a result of their experiences, ‘‘By 1895 most engineers _ 
and architects in Chicago felt themselves masters of their foundation problems. _ 
They accepted the idea that a certain amount of distortion of building — 
due to differential settlement was inevitable, and, if anticipated, was not harmful.” 
Time ultimately has provided the opportunity to judge the relative — 
_ of these early foundation engineers and their treatment of potential settlement. 
Peck cites several cases in which the observed settlements exceeded 12 in. 


_ (300 mm). For the Monadnock Building, the maximum settlement eventually 


exceeded 21 in. (530 mm), and differential settlements in excess of 4 in. (100 


- mm) were observed. Thus, the settlements clearly must have been seriously _ 
underestimated in the original designs. However, Peck reports that aside from 
the three structures with towers, which have been mentioned previously, weal ; 
serious structural damage was reported for 56 of the remaining 57 structures | 
listed in his his study, and 38 were still in service in 1940. Many ‘subsequently 
have been rer replaced by much larger | high ri rise office buildings, but several of - 
_ these early structures, including the Monadnock Building, remain in use today. r 
_ Several basic points are illustrated by the Chicago experiences. First, buildings 
should be expected to settle. Second, differential settlement is more critical — 
than total settlement. Third, structures and their foundations can be designed — 
~minimize differential settlement. Finally, significant settlements can be 


tolerated without i impairing the safety and function of many structures. bs ing 
Deveropment oF Current Crrreria For Dirrerentiat SETTLEMENT 
Almost all current tolerable settlement criteria for buildings are based on = 
the classic studies by Skempton and MacDonald (9) and by Polshin and Tokar _ 
©. Skempton and MacDonald (9) reported observations of settlement and damage - 
for 98 buildings, including both steel and reinforced concrete frame structures 
and structures with load-bearing walls. Buildings supported on spread footings, — 
mats, and piles were included in the study. Damages to panel walls, interior - 
_ partitions, floors, or primary structural members were reported for 40 structures. 
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Because most of the damage appeared to be related to distortional deformations, 


a “concluded that watiie of load- -bearing walls or panel walls in frame structures 4 


= ; is likely when 8// exceeds 1/300 | and that structural damage is probable when _ 4 


8/1 exceeds 1/150. Finally, Skempton and MacDonald suggested 8// = 1/500 


_ The Skempton-MacDonald criteria subsequently were incorporated into recom- 


as a design criterion that provides some factor of safety against any cracking. 


mended limiting values of angular distortion proposed by Bjerrum (1). See Table Sa 
An more recent study by Grant, Christian, and _Vanmarcke G) reviewed - ae. 
wat 
had. ‘suffered some > damage. This “study “supports. ‘the Skempton- MacDonald 
conclusion that cracking should be anticipated when 8// exceeds 1/300. 
_ Both Skempton and MacDonald and Grant, et al. assume that tilting is rigid 
- body rotation which does not contribute to the distortion of the structure, , and, 44 
thus, both have removed the differential settlement due to tilting from the Py 
uted values of angular distortion. Leonards (6) questions the 


to machinery sensitive to settlement 
Danger to frames with diagonals 

Safe limit for no cracking of buildings” 

First cracking of panel walls 
Difficulties with overhead cranes 

Tilting of high rigid buildings becomes visible _ 


Considerable cracking of panel and brick walls 


"Safe limits include a factor of safety. 


this assumption yn for a a framed structure » supported on isolated spread rama 
For this case, he correctly notes that tilting contributes to the stress and “ 
_in the frame unless each individual footing tilts or rotates through the same pe +4 
angle as the overall structure. Because this is unlikely to occur, he suggests 
that the effects of tilt should be included in the differential settlement cr criteria is 
particularly for framed structures on isolated footings. 
pa It is interesting to note that Terzaghi and Peck (12) suggested that “Most 2q 
ordinary structures, such as office buildings, apartment houses or factories can 
th withstand a differential settlement between adjacent columns of 3/4 in.’’ This {a 
is approximately equivalent to an angular distortion of oe if one assumes 
 Polshin and Tokar (8) presented Soviet experiences and included allowable e 
settlement criteria from the 1955 USSR Building Code. These criteria, which a 
int based on more than 25 yr of settlement summarized 


= of and MacDonald. Frame structures a an wells 


3 ‘There are several the Soviet criteria the 
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treated separately. For. frames, the allowable is 
terms of the slope, or the differential settlement between adjacent columns; 


TABLE 2.—Allowable Settlement Criteria: 1955 U.S.S.R. Building 


ype o structure ard clay 


= steel and concrete structures 
_ For end rows of columns with brick cladding since 
For structures where auxiliary strain does not arise during 
nonuniform settlement of 
Tilt of smokestacks, towers, etc. 


Craneways 


atL/H=S wre 


TABLE 3.—Allowable Average Settlement for t Building Types” 


Allowable 
average 
Settlement, 


Bes... 


| Bailing with brick walls, reinforced with reinforced concrete 


i 


which is very. similar to the angular distortion criterion of Skempton and 
- MacDonald without correction for tilt. The limiting values in the Soviet code 


clay 
~ 
00077 
0.003 | 0.003 
| 
0.0003 0.0004 
0.0005 0.0007 
9.0010 0.0010 . 


from m 1/ 500-1 / 200, are same as Mac- 
a Polshin and Tokar’s treatment of load-bearing walls introduced several signifi- 
cant concepts. First, instead of slope, the allowable settlement is defined in 

terms of deflection ratio, A/L. . Then, the maximum allowable deflection ratio 

is assumed to be related to the development. of a critical level of tensile strain 
VE in the wall. For brick walls, the critical tensile strain is assumed to be 0.05%. 

pars these concepts, the deflection ratio at which cracking ooeenean in brick | 


The theoretical relation agrees with observations of damaged and ies 
A buildings. ‘Finally, for /-multistory brick buildings, a _ larger deflection ratio is 


on sand or | r hard clay. Presumably, the : slower rate e of ‘settlement for plastic 
clays allows time for creep of the structure and increases the level of _— 


a 


FIG. 2. —Rectangular Beam for Building (2) 


‘strain, and, therefore, the deflection ratio at which « cracking ais * 4 

In summary, the Polshin and Tokar paper generally supports the Skempton- 
MacDonald criteria for framed structures but imposes more stringent limits on 
differential settlement of load-bearing brick walls. Polshin and 1 Tokar also — 
introduce the important roles of critical tensile strain, | (L/H) as ay measure 

of structural flexibility and the rate of settlement. The concepts of Skempton 

and MacDonald and Polshin and Tokar can be found in most current settlement 


The state-of-the-art report by Burland and Wroth (2), , presented at the 1974 
_ Conference on Settlement of Structures at Cambridge University, has contributed | 
= to the understanding of factors that influence the allowable settlement 
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7 7 value of tensile strain is reached, Burland and Wroth proposed that a building | 
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may be represented by a (Fig. 2). This” 
model then is used to show the factors that affect the settlement at This 
7 
In this idealization, the deflection ratio, A/L, at which the | critical tensile 
‘strain, €.,ix, iS reached at some point in the beam is used as the criterion for — >? 
_ allowable deflection or settlement. The critical tensile strain may develop either 
by direct tension in the extreme fiber, _A, or by diagonal tension along the 
‘neutral axis, B or B’, (Fig. 3). The direct tensile strain, cs and the diagonal = 
; - strain, €,, are defined in Fig. 3 in terms of the bending and shear stresses, 
respectively. Using elastic theory for deflection of beams, the deflection ratio, 
4 / L, and the slope, or angular distortion, of the deflected beam can be expressed 
in terms of the elastic properties of the beam and either of the critical = 
strains, €, or i Burland and Wroth considered the deflections of a simple 7 ) 
beam, incheding both bending and shear effects, for | "two loading conditions: 
@ A point load, P, at the mid-span; and (b) a uniformly distributed load, w. — 
® relationships between A/L and critical tensile strain for these cases are a 
‘The role of various factors on the establishment of the allowable (A/L) may 
I be observed by plotting the relation of WE to L/H. Fig. 4 shows these _ 
relations for the case of an isotropic beam with the neutral axis at the mid-depth 7 = 
. and Poisson’s ratio equal to 0.3, for which the ratio of Young’s modulus to _ 
> shear modulus, E/G = 2.6. Curves are plotted for the central point load - 
and the uniformly distributed load assuming either bending strain or diagonal 
strain may be critical. When bending strain is assumed critical, the curves are 
relatively insensitive | to the assumed load distribution. When | diagonal strain 
_ is assumed critical, the curves for the two loading conditions are again similar 
if shear stress at the quarter- point of the span is considered. However, if shear 
_ at the end of the beam is considered, the results are quite different for the 
two loading conditions. Burland and Wroth suggest that the theoretical shear — 
stress conditions at the quarter-point are more representative of practical 
conditions in real beams. Then the load distribution is relatively unimportant, — 
and the concepts for leniting Values ¢ of (A/L) may be developed using only 
the case of the central pointload. 
_ Fig. 4 clearly shows the effects of the critical tensile strain, €.,,, and the — 
_iL/#H ratio of the beam. For a given value of L/H, the limiting A/L will 
be directly proportional to ¢€.,,,. For low values of L/H, diagonal strain is 
critical, and for large values of L / H, direct bending strain is critical. The allowable _ 
4/L increases as L/H increases and the beam becomes more flexible. ea 
Burland and Wroth also suggest that the effect of a mat foundation, which — 
is very stiff in the lateral direction, might be / modeled by assuming the neutral | 
axis is at the bottom of the beam (Fig. 5). For the usual sagging concave 
upward) pattern of settlement, the entire beam will be in compression, and _ 
Pt the diagonal strain condition will be critical for all values of L/H. However, — 


if the settlement curve is concave downward (hogging), i.e., the settlement 
of the ends is greater than the settlement of the center, the direct bending 
strains become critical, except for small values of L/H. For this condition, 
the allowable deformations will be significantly smaller than for the case with 

- the neutral axis at the middle. This result suggests that the allowable differential 

settlement should be smaller for cases in which the settlement curve is concave 


q 


cuadiieeat a e.g. in the case of subsidence of ends due to o adjacent excavations 
or differential heave due to swelling soils. bas saves! 

' _ Burland and Wroth also note that many structures are more flexible in shear 
than in direct tension. They suggest that this condition may be modeled by — 

_ varying the ratio E/G of the rectangular beam. The results show that diagonal 

Strain _ becomes more critical as th the beam i is assumed more flexible in shear 


cases, illustrated in in Figs. 6; and 7: 
|. Diagonal strain will be oii: for framed : structures, which ‘iii are 
ty: res flexible in shear, and reinforced load-bearing walls, which are relatively _ 
_ stiff in direct tension. This case may be approximated by the behavior of the 


: ¥ rectangular beam with a al equivalent E/G, a. shown woe curve 3 in = 
6 and | curve | in Fig. 


(hogging only) 
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5.—Relation of to for | Axis at 


= strain will be critical for masonry walls and structures, 
which have relatively low tensile resistance. For the sagging mode (settlement 
curve concave upward), the behavior may be modeled by the isotropic rectangular 
- beam with the neutral axis at the mid-depth, as shown by curve 2 in Figs. | 
- 6 and 7. For the hogging mode (settlement curve concave downward), the behavior — 
may be approximated by the rectangular beam with low equivalent, E/G, and 
the neutral axis at the bottom, as neennes by curve 4 in Fig. 6 and curve 


oth 


_ The curves Soy A Lane ip in Fig. 6 are converted to curves of A/L in a Fig. 


crit 


cracking of eareinforced load- -bearing walls, which assumes = = 05%, 


tho 


© allowavle L huch larger tnan lor une Cases With tne lower 
values of E/G, and A/L increases rapidly with increasing L/H. 
o er... the basis of this simple beam analogy, Burland and Wroth propose that 
: limiting deflection ratio criteria should be developed for at least three different | 
| 
‘ 
Bottom 
a 
7 
| 
and the Skempton-MacDonaid criteria 0 UU an 
_ In order to relate angular distortion to deflection_rat 


assume that the maximum angular distortion of the beam ‘the 


_ The curves in Fig. 7 are - compared with observed settlements and damage 
of real structures in Fig. 8. For frame buildings with L/H less than 3, the 
E/G = 0.5 very stiff in shear 
E/G = 0.5 neutral axis at bettom 
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6. —influence of Equivalent on of 


--neutral axis at mid-depth, 
bending strain critical 
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FIG. 7. Criteria for Limiting Values of 


provide and satisfactory limits for ‘the “deflection ratio. For L/H 


greater than 3, the Skempton-MacDonald criterion is the more. coaservelive 
ror 
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of the two. However, no field observations were available for comparison. 


m ; For load-bearing walls in the sagging mode [Fig. 8(b)], the Skempton-MacDonald © 


criterion appears to be unconservative while Polshin and Tokar’s limit and curve ~ 
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2 to ‘reasonable ‘limits. Finally, Fig. 8(c) shows that 


- -Joad- -bearing walls in the hogging mode may crack at deflection ratios that are 
4 much smaller than the Polshin-Tokar limits. However, curve 3 appears to provide - 


| 
q 
(c) Hogging of load bearing walls fp 
FIG. 8.—Comparisons of Obser 


reasonable limit for this case, but relatively few data are available. 
4 In summary, the Burland-Wroth beam analogy provides a conceptual mecha- 
nism for understanding the factors that may influence tolerable settlements. 
The model can be used to demonstrate the. effects of the « critical tensile strain 
_of the structural materials, the L/H ratio of the structure, the relative stiffnesses 7 
_ of the structure in shear and direct tension, the longitudinal stiffness of the 
foundation and the mode of the settlement curve. It also provides some insight 
into the type of information that that eel be included in meaningful settlement - 


Summary a1 AND Concuusions: 


“of ‘architectural elements. Angular distortions of 1/ 150 c can n be peta to cause 
"structural damage. For load-bearing walls, the tolerable differential settlement 
= smaller than the values for framed structures and is influenced by the L/H 
ratio of the wall. The allowable differential settlement i is smaller than the values 
i framed structures and is influenced by the L/H ratio of the wall. The 
- allowable differential settlement also is significantly y smaller when the settlement 
_ Future improvements in establishing tolerable "settlements for various types 
of structures will depend on the continued accumulation of reliable settlement 
records for both damaged and undamaged structures. The beam analogy of 
“contd and Wroth (2) provides a conceptual understanding of the factors that 


control tolerable settlement. These concepts can be used to identify information ; 
required for meaningful interpretation of case studies and to separate data an ; 
case studies into appropriate categories of structures. __ 


ApPENDIX —Dertection Equations ror RECTANGULAR Beams 


i 
‘The deflection ofa semnnaitis beam due to bending and shear can be computed > 


“from elastic theory. The direct tensile strain, €,, in the extreme fiber and the | 


diagonal tensile strain, at ‘the neutral axis is also can computed from | elastic 


_ For a simple rectangular beam of ‘span L, cross-section » #xi, él with 


the neutral axis at the mid- id-depth oe 
| 


n = = diagonal tensile tensile ‘Strain at the neutral alaxis . 

in which M = bending moment; V = shear force; H = depth of beam; L = length — : 
of beam; E = Young’s modulus; G = = shear A =H x | = cross-sec- =Hxl= 
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| 

_ The current criteria for tolerable settlement of structures are based primarily 

on the works of Skempton and MacDonald (9) and Polshin and Tokar (8). Different | 
4 _¢riteria are required for different types of structures. For frame buildings with bod 
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c=H/2= distance from neutral axis to extreme fiber. 
A relation between deflection ratio, A/L, and angular distortion 8/7 = 
fi i may be established by assuming that B can be represented by the slope a 
some portion of the beam. Three definitions of B are considered: Bo= = i | 
at the support = maximum slope; ,. ta = Slope at L/4, the q -point of the 
span; and Bavs = = slope of chord between the support and wit 


and 

The corresponding shear and d diagonal strain are ai novig 2 baw 
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in 1 which - and C, are coefficients that depend ‘on the definition of B. Vanes 
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Ld | The following : symbols are used i in this paper: 


A = _CTOss- -sectional area of a Tectangular beam; 


| — 
Whos 
§ 
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shear modulus; 


depth of rectangular bens: 
= moment of inertia about neutral ax axis; “2 


Young’ s 


L = span length of 
= bending moment; 

vertical shear; 

uniformly loads 

= relative deflection; 


shear strain; 


vertical 

= tensile stress; 
shear stress; and 

tilt or rigid body rotation. 


Subscripts 

crit = critical value; ke 
= point along span length; 


ath along span length. 
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LAGOVEN’s -Amuay Refinery lies on the west coast t of de 
_Paraguanda by the Golfo de Venezuela. The peninsula lies between \arncas and 
Maracaibo, about 200 km east of Maracaibo, 
Along the Amuay coast exists a cliffside having a typical height of iiss m- 1-25 
m. LAGOVEN (and its predecessor, Creole Petroleum Corp.) constructed earth — 


embankments to enclose indentations in the cliffside, termed quebradas, for : 


the storage of fuel oil. The Amuay Refinery has three of these oil storage 
reservoirs, FORS-1, FORS-2, and FORS-3. Lambe (1956 and 1963) has aaeenieadl 7 
these open reservoirs for storing fueloil, 
_ From industrial and domestic activities, perched water has developed ai 
a stratum of weak clay, , brown fat clay. The development of this perched water 
has triggered landslides’ along. ‘the Amuay cliffside. Much of the cliffside at 
-Amuay | exists in an unstable condition. Landslides have occurred at a dozen 
_ locations; surface tension cracks indicate incipient slides at many other locations. 
Landslide means ‘‘downward motion of a natural slope due to gravity.”” This — 
_ paper uses the term /andslide to describe a shear slide along the cliff even 
_ though water aids gravity in actuating slope movement and some of the Amuay 
cliffside has received alterations, smoothing, slope protection, etc. 
_ The consequences of a landslide along the Amuay cliffside can and do v vary 
F considerably. A slide in a relatively remote location merely causes inconvenience 
a the effort required to dress the failed cliffside. At the other extreme, + 
_ Slide near a critical refinery facility could have catastrophic consequences. For 
_ example, a landslide which breached an oil storage reservoir could spill ae 
~ than 10,000,000 barrels of oil into the sea, ‘resulting i in a loss of several hundred © 
_ millions of dollars worth of stored oil, monumental political problems both national — 
_and international, and large environmental damage. 


_ Fornearly 30 yr LAGOVEN (South America) has supported and helped execute © 
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particular, a landslide near an important facility. During this time, we have 

performed tests, both in the laboratory and in the field; made stability and 

flow analyses; taken field measurements (especially measurements of pore 
g Pressure and movement); and carried out preventive ai and remedial project s. 
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publication describes the geotechnical safety pr program. _ This prog program rests on 

Predicting instability and then taking appropriate -remedial measures. 
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_AMUAY CLIFFSIDES 1507 
very well. Essentially no problems have developed with 
man-made embankments of compacted soil. Difficulties, however, have occurred 
at ‘many locations along the cliffside, including | the natural earth slopes enclosed 
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in the reservoirs. In fact, oil stored in a reservoir applies a a force which helps 
resist a slide of the inside slope of the reservoir, = = emossty 


Fig. 1 also chronicles geotechnical malfunctions, measures taken to prevent 
landslides and measures taken to repair lastdelides. Fig. 1 indicates that the 


= in on storage reservoirs. Stadility prodiems exist only along the 
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Th This section examines the mechanism m of a landslide at Amuay. The following - + 
three factors have overriding importance: (1) Geometry of the moving land 


This ‘section takes an initial look at these factors and the following sections 


examine eachone more closely, j= 

of the landslides i in the Amuay Refinery have involved the movement 

fal one imagines a wall between two wedges, the wedge to the left of the wall 


_ corresponds to an active wedge and the wedge to the right of the wall _ 


Fig. 2(b) shows a freebody of the passive wedge and Fig. 2(c) lists the forces 
- on the passive wedge. As the figure shows, the actuating forces tend 


resisting forces tend to prevent | the | 
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outward movement. In analyzing a ‘Stability : situation, , the engineer determines 


a landslide the force T equals | the maximum ‘resisting force that the soil strength 
can generate. Factor of safety equals the ratio of maximum resisting force 
toactuating shearforce. 
Fig 3 shows the soil profile at Amuay and presents the heads existing in _ 
ground water. The Stratum of fat clay serves as a seal t to entrap 

water entering 
 piezometers inserted at various locations. in the profile. Fig. 3(b) portrays the 
heads existing in the ground water. As Fig. 3(b) shows, the development of 
perched water increases the pore pressure on the top of the brown fat clay. 
; _ As later sections in this paper discuss, the buildup of pore ee on = 
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An early step in analyzing a ‘slide consists of « detaining the geometry of 
the earth mass which slid. Using visual observations and common —— 
the engineer generally determines the exposed boundari 
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gim-2) 5.—FORS-3 Slide, Dec., 1976 


_ Aided by borings, excavations, surveys, probes, visual examinations, soil tests, Z 


etc., we have tried to define the geometry of 10 slides at Amuay. Figs. 4, 7 
- 5, and 6 show the results of six of these attempts. Fig. 7 shows an actual» 
_ slide at Amuay. On these geometries we have: wall 


 ¥ 1. High confidence in: vertical crack; breakout of slides; and face of slope, 
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— confidence in: teesdion of the shear surface in clay; and angle 
t 3. Little confidence in: location of the surface between the active and passive 
On two landslides we have located with good accuracy the shear surface 


e, 1979 


in a very thin i mm thick) 

rone exi: baler aan m below ‘the top of the fat clay. On another slide we 
tt cated two thin shear surfaces spaced about 300 mm apart in the stratum 

_ We have not located on actual slides any surface acting like a retaining wall 
a between the active and passive wedges. Further, we lack field data which would 
_ justify one to postulate some sort of oe movement hove: ve failure 
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failure geometry, not circular f failure geometry, for subsoil conditions 
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FIG. 8. —Geotechnical History of Amuay Cli Clift hes 


evaporation of ground w water r during recent geological ti times, caps the Amuay 
, it cracks before © 
landslides 0 occur ir and furnishes little or no resistance to : an wn incipient landslide. i 


A 2-m- 5. -m stratum of caliche, consisting of precipitants produced by the 
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_ water during the Tertiary Period. Extensive and well-preserved malacological 
fauna indicate deposition during the Medium and Upper Miocene and PI: ‘ocene _ 
Epochs (Gonzalez de Juana, 1938; Feo- -Codecido, 1971). Most 
apparently came from the mountain range directly to the south of Amuay 
=> ~Codecido, 1971). The maximum overburden on the present surface sediments 
i mill during the Pliocene Epoch; geologists have estimated that 100-200 ! 


m of additional sediments existed at that time. Results of oedometer tests suggest 
that about 80 m of additional sediments existed above the present- day ground — 
surface at the refinery site. During | periods when a low sea level existed, erosion — 


_ of these sediments occurred. Desiccation of the remaining sediments could have - 
extended 20 m or more below the present ground surface. As a result of erosion | 
- and desiccation, the sediments at Amuay have experienced overconsolidation, 
: probably in excess of 1,470 kPa. The soil below the caliche varies —s 
in both iateral and vertical directions, es from sand with lenses of clay 
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to silty sand to silty clay. We term this variable layer as silty sand. In the 
typical landslide along the Amuay cliff, about one third of the failure surface _ 
goes through the silty sand. Below the silty sand lies a stratum of brown fat — 
clay. About two thirds of the shear surface of the typical landslide lies in 


principles of the Stress Path Method (Lambe and Marr, 1979). As best we 
could, we traced out the stress history of three soil elements of interest on 
a potential landslide failure surface. Fig. 8 shows the results of this exercise. 
' We studied eight events in the history of the soil elements and identified stages 
at the end of each event. For each of these ‘Stages we estimated stress values 
for the three elements under consideration. Fig. 8 tabulates the estimated stresses. 
Fig. 9 portrays the stress paths for element B, which approximates the average 
_ Having the stress paths, we then conducted two series of tests. For one 
Series we sedimented fat clay from a We then n subjected 
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ve the s sedimented clay to the stresses given in Figure 8 and sheared samples a 


In the other series we obtained undisturbed samples of fat clay 
_ from various locations in the refinery. From one location we obtained a sample 
* at Stage C, from another a sample at Stage P,, and another sa sample from Stage 
. We ran drained tests and undrained tests on these samples. i en 
Zz “We obtained values of residual strength two ways—one set of tests on an : 
q ‘undisturbed sample from the shear surface of an actual landslide and + ee 
of tests on samples subjected to large deformations, 
Fig. 10 summarizes the results of laboratory tests on the undisturbed samples _ 4 
and the s samples taken to residual strength. The residual moh depends upon = 


of the clay and, further, that the test resus don not give straight- -line 


Perched Water (intact slope) 
Perched Water (cracked slope) 


. 10. —Strength of Brown Fat Clay 


loss contributes to described by Skempton (1964, 
Pore Water | = 
As Fig. 3 water exists above. the brown fat clay. ‘The perched 
Y _ Water causes positive pore pressures in the upper portion of the ‘fat clay and 
silty sand overlying the clay. Fig. 3 indicates a pore pressure in the fat clay 
at elevation +12 m equal to 98 kPa. Without perched water, the water in the 


fat clay would exist at a A negative pressure. if 


= the clay. These contours indicate that the perched water layer exists over most, 
if not all, of the refinery area. The contours further suggest that the source 
of water lies “a the len at Field measurements dating back to the mid- 60s 
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that total head | has during: the last decode... 
a The head described to this point comes from piezometers having sensors ; 
either above or below the fat clay. The heads represent values associated with 
seepage in the soil above and below the fat clay. Fig. 12 suggests the distribution 
of pore pressure in the fat clay near the southeast wall of FORS-3. The pore 
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. FIG. | 12. 2.—Pore ‘Pressure i in Brown Fat Clay, Southeast Wall, FORS- 3 ose 
"pressures plotted in Fig. 12 came from measurements made on a pore pressure 
We have not yet explained some of the pore pressures shown in Fig. 12, 
aa the very high pore pressure (207.9 kPa, giving a total head of +30.2 
-m) occurring at elevation +9 in the fat clay. In making slide analyses and — 
Stability predictions, we ha have used used values of of pore to 
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existing above ‘the fat — reason this pore exists 
_ near the surface of the fat clay at the start of a landslide. What little evidence — 


s we have suggests that the shear zone in the fat clay has such a small thickness" 


that drainage probably occurs during the early stages of a landslide. We ‘suspect 
that the high pore pressures measured well within the fat clay result largely 
_ from excess pore pressures developed by deformations occurring during a 
landslide. We have low confidence in the correctness of this explanation. During | : 
_ the years ahead we shall focus attention on pore pressures within the fat clay = 
in hopes of establishing the magnitude of pore pen: which hich will exist at 


the start of a landslide. 


_ Sliding occurs when the sum of the forces eG Dab becomes equal © 
to the sum of the maximum forces reg movement. a to Fig. 2, 
q TABLE 1.—Stresses for. Element 
pascals | s | line 


17.7 1,373.0 


Increase P’ (placement of surcharge). 


Increase U, ee -up of perched water, development of excess pore pressure 


. Increase U, (buildup of ‘perched wate water, development of excess pore pressure — 


Decrease c’ or (deformation of the clay, thereby lowering the strength 


6. Alter W so as to increase G or decrease W-U,,orboth, 

_ In the early 1970s we prepared Fig. 13 from analyses of actual slides and _ 

from saaeeed tests on undisturbed — of soil from —_ shear zone of 


_ _we see that any one, or some combination, of the following can initiate a slide: _ 
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e east borrow slide —_— We tried to relate | o’ with age « ge of slope. We re reasoned» 


th 

that a decrease of ’ with age triggered landslides. This approach did not work 
1. We had trouble determining the age of a slope. itgid 
2. Time per se probably does not constitute the primary variable. got! 


vil One cannot afford the flat slopes obtained from a design using a value = 
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FIG. 14.—Strength and Stress Paths for ElementB 
va Using the strength lines and stress paths for Element B shown in Fig. 14, 
we can show conceptually the factors causing a slide and can indicate the relative a 
_ importance of each factor. Fig. 14 gives strength and stress paths for Element — 
B, which approximates the average element for the shear surface of the passive 
_ wedge. Table | tabulates stresses for Element B for Stages B, P,, and P. 2 
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can see that going from — P. , for example, involves three 
Gaper oF » case invelv ot efiffeide 
= vw An increase in n applied shear stress (q goes from +117. to kPa——83.4 kPa, 


ee Y decrease in effective stress (p’ goes from +1,373 kPa tc to +157 kPa, 


3. A decrease in strength line (for any p’ value, 4 loses about 3/4 in going 
from Stage C to and and drops about 1/3 i in /3 im going C to Stage 


The rise in perched water causes all of the increase in applied shear stress, a 


alll of the decrease in effective stress, and a major part of the drop in strength 
7 from Stage B to Stage P.. The strength drop from P. tol R results: from deformations 
caused partly by the i increase in perched water and partly by « other factors. 
"Deformation in the clay, e.g., from the cyclic filling—emptying of the oil in a 
reservoir contributes to a decrease in strength. ‘perched water 
tas Excess pore water pressures which develop from the undrained or partly 
‘drained, loading or unloading of the clay can contribute to the initiation of 
‘the landslide. _As noted earlier in this ‘Paper, have not sorted out 


headings 


culate 


We can sinha the ti tasks facing the geotechnical e engineer solving stability 


problems into one of two categories, namely: 


The slide analysis, a slide has occurred, L.e., the factor. of safety 
The Shear stress along the failure surface the shear strength of of the 


Ey. geotechnical engineer does a slide analysis to gain insight and an 


use a 1 slide analysis to \ validate a prediction method. _A slide ° analysis constitutes — 
an autopsy—i.e., a Type C Prediction according to Lambe (1973), is 
- stability prediction constitutes an assessment of stability—expressed 7 
a factor of safety or the probability of failure—for an unfailed slope and 
given set of conditions. This corresponds to a ‘Type A ‘Prediction. One may 
any one of the following four situations: = 


- 
summary, the development and increase of pe 


18 
4. Future slope, future conditions. a 


Ina a stability Predict i », te engineer 1 must forecast what events can initiate 
sliding. and usually he st 
exist at the future condition. Stability prediction usually presents the engineer 
with more complications and uncertainties than does slide analysis. This paper — 
“Fig. shows an analysis of the FORS- 3 Dec., ‘1976 slide. wot 
We determined the geometry of the sliding mass as discussed in the section 
on geometry. We took for the ‘silty sand equal to’ 33° based on ‘triaxial 
and direct shear tests results (the residual strength approximately ‘equals peal 
peak strength for the silty sand around FORS-3). We selected the strength 
of the clay from Fig. 10 using Stage P. for the original landslide and residual 
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15. of Dec., 1976 Slide, FORS- 


Fig. 1s the of ‘total head on of the fat clay used in 
_ The analyses shown in Fig. 15 prove that—knowing ahead the answer—we 
can calculate the correct factor of safety (unity) for a slide. This fact indicates _ 
either we used a correct procedure and correct parameters < or a fortuitous 7 
cancellation of errors occurred in the analysis. We suspect the latter since we ; 
made several simplifications in the slide analysis ce. Ee 
» bos aa Yo a vis 
' Using pore pressures from a steady-state flow net but measurements a . 
a slide show much larger values. 
2. Assuming the inclination of R,, 
3. Assuming deformations occur along two rigid wedges whereas the measured 


ia _ position of surface of slope Suggest a more complex pattern of deformation. _ 
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; measures, the instability conditions of the cliffside will continue to deteriorate. — 
_ The consequences of a landslide along the cliff depend upon the feclins 
4 near the landslide. A slide which breaches an oil storage reservoir could dump 


4 millions of barrels of fuel oil into the Caribbean Sea, an undesirable consequence 


ra The instability exists in the natural cliffside. No stability problems have arisen 


in the man-made dams which enclose quebradas to form the oil storage reservoirs. © 
The failure surfaces of the landslides lie mostly within a layer of weak plastic 
clay, termed fat clay. Most of the landslides involve wedge-type failure geometry. 
_ This paper presents test data to show ‘that the strength parameters of the — 
Perched water arising from | activities within the r refinery aye developed above = 
the stratum of fat clay which acts as a seal. This perched water initiates the — 


~ landslides. The rise in pore pressure from the the poses) water causes: 


le ey in effective mavens, ant thus soil 


¢ A drop in strength parameters, arising from deformations : ah 
3. A rise in shear stress along a potential landslide shear surface onal 


wedge- “type s slide analysis pore pressures inferred from readings 
tests 
“4 yields a 4 factor of safety of | unity. The pore pressures used correspond t to , pressures 
existing prior to a landslide. The shear strength selected depends upon the 
Stage of the clay just prior to the landslide | 
The slide analysis illustrates that—knowing ahead the answer—we can calculate 
the correct factor of safety (unity) for a slide. This fact indicates either we 
used a correct procedure with correct parameters or a fortuitous cancellation | 
pa paper - employs laboratory and field measurements a 25- year 
"period. Many people contributed to the collection of these data. The writers 
have received support and encouragement from LAGOVEN engineers, especially — 
a Carlos Hernandez and Mario aren. The writers S appreciate the contributions _ 
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By ‘Leland M. Kraft, 'M. ASCE, A. Focht, F. ASCE, 


“average frictional resistance, engrossed as a fraction of the average ope 
shear strength or average effective overburden pressure, decreases with increasing 


that are the basis 
the studies include: “() Tests with variations in set- “up time and stress history; 
_ (2) tests on piles redriven to deeper penetrations after earlier testing; and - 
tests on piles with oversized closure plates. These factors can combine to amplify : 
the influence of pile length on axial capacity, q 
_ The influence of the length effect suggested by earlier studies has significant 
economic impact on offshore structures that are supported by large-diameter, 
open- -end pipe piles driven to penetrations of 200 ft-500 ft (61 m-152 m). 
7 Consequently, controversy has continued about the applicability and validity 
of current capacity computation methods. In the past two or three years a 
- modest amount of new load test data has become available, and several model 
7 - and analytic studies have been completed. Therefore, to place the length effect 
into proper perspective, the writers reexamined the influence of ’ pile length 
on frictional capacity in clays using current technology and philosophy. Some 
3 of the methods now being used to compute pile capacity and to study ol 


that affect pile capacity are summarized first to set the stage for this —" a 


Methods used to evaluate and compute unit ‘shaft friction, f, for offshore 
_ piles in clay include the Lambda method, (32) the methods recommended by — 


API RP 2A, (1) the B methods, (6,15,24,27) and the Janbu (19) method. The = 
: API methods may be referred to as the f = s, (s, = undrained shear strength) 


4 Mgr., Special Projects Group, McClelland Engrs., Inc., Houston, Tex. 
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approach for highly plastic clays or the a approach for stiff, overconsolidated 
clays of low plasticity. The B method and Janbu’s method are effective stress — - 
4 approaches. Other effective stress “approaches are ‘being developed, but 
‘not yet ready for design applications (14,21). The Lambda method, API methods, 
_ B methods, and Janbu’s method are summarized in Table 1. A length effect 


TABLE 1.-—Design Methods | 


 &,,, = mean G,, along pile Ab 


RP 2A (1980) | f = B, = 


Oc:°a=1, but f+ the larger of | or (s, 


(ii) to plasticity clay te sone 


1-05 din. var")* a 


0.5 


(1973) NC: B =(1- sin) tan > 

= effective angle of internal friction _ 
Meyerhof (1976) | NC: B = Func’ 
B= 1S sind) tand VOCR 


ijayvergiya and 
_Focht (1 (l 


which p = tan ®/|r|;7 = roughness = = Func” 
(L); and a = attraction 


Func" (x) = Function of x; lin. le 


is ‘implicitly included in the Lambda, Flaate and Selnes, and Janbu ‘methods; 
_ Meyerhof incorporated a length effect into his chart for normally or slightly FA 
clays; and application of APL Method 2 produces a “length 


effect for normally consolidated clays. __ 
of Methods. Clays: Predicted B profiles 
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DRIVEN RIVEN PILE FRICTION 
i< (eatio of average f to initial effective overburden pressure) for several of the — 
pile capacity methods are compared in Fig. | to show the implied length eft fects. 
B values from the ‘method are obtained from f = + 28 in which. 
B= the ra ratio of undrained shear strength to effective overtesden pressure 
_ (s,,/F,.); F ym = the mean effective overburden pressure over the pile penetra- _ 7 
= tion; and \ = an empirical coefficient. Because the \ method was not intended 
for short piles in normally consolidated clays (32) and overpredicts capacity g 
for less ‘ft, the band for is not plotted wil 
shy poo : 0.2. 02 04 05 


= 
io 
= 


FATE 
AND 


NORMALLY 
CONSOLIDATED 
—Comparison of Values (1 ft = 305 m; 1 ksf = 47. .9 KN/m* ) 


average B value is as low as 0.1 for a B, value of 0. 
- Measured and computed capacities were compared Statistically for the Lambda 

API methods by Drewry, ‘Weidler, and Hwong (12). The means of the 
ratio of calculated to nian eens were — the A method and 1.29 


differ by about +25% relative to the median of the four methods. 
@ The average B values implied by the API methods are not included in Fig. | 
_ 1, but in normally consolidated clays, B ranges between 0.2-0.4 for all penetrations 
r for Method | and for pile penetrations of less about 60 ft for Method 2. For 


___Overconsolidated Clays. —Drewry, Weidler, and Hwong (12) found that the 
- ‘fatio of calculated to measured capacity in overconsolidated soils was 0.83 for 
the API methods and 1.04 for the Lambda method. The standard deviations 
were nearly the same Ld both methods (0.18 ; the Lambda method and 0. 20 7 
for the API methods). > 
a Pile capacity methods. fe or overconsolidated clays ys are more. difficult to compare = 
in a generalized way than methods for normally consolidated conditions. There- 
 fore,a comparison is made for a specific site condition. Computed pile capacity — 
curves for an _overconsolidated soil deposit in the North Sea are cmeued 
in in Fig. 2 for four methods. The deviation the average of the curves is 


20 


z 
2 
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FEET 


FOR JANBU APPROACH 23 DEGREES, a= 1.29 KSF 
OEPTH AND 0 AT GREATER DEPTHS 


= 4.45 KN: = 47.9kN/m*) 

“about +60% at 50 ft (15 m) penetration and shoe at 250 ft (16 m) 


Penetration. In terms of required pile length, the difference i in the pile penetration y 


Primary Causes.— Although a length effect does exist, it may be less pronounced 
than indicated by Fig. 1. Following are the four primary factors that are thought s 


alterations to the soil fabric, and hence stress-strain strength 
properties, with increasing pile length passing a soilelement. = = 

; ” A strain or displacement softening response of the soil during pile loading. 


After reaching a peak shear emote, te stress transfer decreases with additional : 
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the API methods, and the standard deviations were 0.31 for the method 
_and 0.44 for the API methods. This suggests that the +25% deviation between 
| 
q 
| 


D VEN PILE FRICTION 


‘The in influence of thr of three- dimensional (e.g. surface conditions 
effects, Jack of plane strain | or plane : ‘stress conditions pile 


Once several feet of pile pass a point, further alteration to “— soil fabric 
3 = to additional pile passing a point should be minimal. Additional changes _ 
= in the soil fabric will occur during consolidation after « driving which : should — 
have a consistent pattern along the pile length. After the consolidation, soil 
= fabric along the pile will thus probably be relatively consistent if the soils were 
originally similar along the pile; 
Three-dimensional considerations at the ground surface and at the pile tip 
may not be important to the capacity | of long piles. Another three-dimensional 
effect is lateral pile movements during installation. Lateral pile movements can 
- create an oversize hole which forms gaps between the pile and soil when piles 
are driven into strong, overconsolidated soils (31). If gaps develop, the load 
_ transfer is reduced. Even if a gap is not developed, the lateral movement can — 
- result ir in low load transfer for piles in stiff clays because only low lateral ‘total . 


appear to diminish with increasing pile length. 
_ Aurora, Peterson, and O’Neill (2) report capacities of 1 i in. (25. 4 sta siiiniae 
_ piles driven 38-100 in. (965-2,540 mm) into an overconsolidated compacted clay. 
These results indicate that the ratio of load transfer to shear strength is less — 
in the upper 30-50% of the pile length _ compared with deeper sections. From © 
our interpretation of other data in the literature, the effects are less severe _ 
in weak deposits (31) and may be less severe when piles are pushed rather 
than driven if alignment is maintained (7,29), 
Variations in soil properties along the pile may also be important to the maximum 
shear stress that the soil can transfer to the pile. Piles driven into a profile 
of sand over clay can drag down “sand between the pile and clay to provide 
‘shaft friction greater than that of the clay alone (31). Seasonal changes can 
result in shrinkage and swelling to further complicate interpretation of pile capacity 
Ponce especially for short piles in heavily overconsolidated a 
i Undoubtedly, the four primary factors that contribute to a length effect are 
interrelated, and all affect pile ‘capacity, but the major contributors are thought — 


s be the strain (displacement) softening phenomenon and lateral pile movements 


Soil Softening. —The load transfer, as defined by shear stress-pile displacement 


_ (t-z) response, may exhibit a displacement softening behavior as shown in Fig. | ; 

| 3. If displacement softening occurs, the peak pile capacity corresponds to the 
capacity associated with the peak load t transfer only when the pile is rey 
stiff and in a homogeneous soil. Short piles tend to have greater relative axial 


is likely to be more pronounced in long piles than short piles. It may also 
_ be more pronounced for large diameter than small diameter piles for reasons | 
presented subsequently. Fig. 4 shows the relative resistance that is developed _ 
at ultimate load for a less than infinitely stiff ple. tt 
_ The effect of displacement softening on ultimate pile capacity has been <_< 


stiffnesses than do long piles. Thus, the displacement softening ZS 
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recognized, and t-z _z analyses can be used to the effect for a 
¥ set of i input ¢ conditions. The pile and soil properties that influence the development 
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FIG. 4.—Load Transfer Along Pile 


of ultimate pile resistance can be illustrated with an idealized t-z curve. Murff 
(26) used tl the curve in in Fig. 3 with ith = to compute the of relative 
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PILE FRICTION 


pile- -soil on the ultimate pile capacity ofa pile ina 
soil. The ratio of the ultimate capacity to the capacity if the peak transfer 
- is mobilized simultaneously ateng the pile, is related to the amount of softening, _ 
and the pile-soil stiffness, 7, , defined by Df,,,, L’/(AEu*), in which D = the 
pile ee; L= the pile length; A = the cross- sectional ar area of the pile; E= 
max = the peak soil-pile friction; and u* 
3 movement at which f,,,, is developed. Murff’s results ; are e plotted in Fig. 5 
as three solid line curves, in terms of a mobilized friction ratio f,, /f,,,, and — 
pile-soil stiffness ratio for three valuesof€. 
Other einen? were performed as part of this study to evaluate the influence 
MOBILIZED FRICTION RATIO, - 
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5. —Mobilized Pile Capacity with Changes in Pile-Soil Stiffness 
—_ of the distribution of | the peak : shaft friction, f,,,,, along the pile, the degree a 
of strain softening, and the variation of u* with depth. These supplemental _ 


results are included in Fig. 5 as individual points and a dashed curve. The — 


q 


and u * at the 


max 


“soil stiffness Tatio is defined in terms the average 


capacity to the capacity ‘if peak transfer is s mobilized is influenced 
_ primarily by (1) Pile-soil stiffness; (2) amount of softening, €; and (3) rate of 


_ softening, B. — results clearly demonstrate that it is unlikely that the sine 
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/ 
a in proportion to the reciprocal of the depth. The influence of the distribution Re 
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friction for a pile can be reliably related to a single para parameter « obtained by. 
combining all or eliminating some of these factors. 
Critical —Four key parameters &,u*,p) are rey to define 


- load transfer, can be determined with the sonata described by Kraft, Ray, 
and Kagawa (22). Effective stress concepts offer potential to predetermine the 
 enduaiaa: of maximum load transfer, f,,,,, and the influence of lateral pile 

_ movements during installation on this transfer (21). The amount of softening, : 
7  €,and the rate of softening, » , are not a readily amenable to theoretical predictions. a 
Some insight into the magnitude of € and p can be gleaned from measured 


t-z responses and from laboratory tests designed to simulate load transfer —— 


The soil next to a pile is severely remolded during Siihetiation, and the driving 
_ process followed by the consolidation process may result in preferred particle 
alignment parallel to the pile axis. The magnitude of strain softening of this 
remolded, reconsolidated ‘material probably differs from that of the 
material. Typical values of gE from measured t-z responses and laboratory 
simulation studies have ranged from about 0.8-1.0, although smaller values may 
be encountered in some cases (8,11,23). Factors such as overconsolidation ratio, 7 
_ soil sensitivity, relative density, and grain size characteristics may affect ‘oul 
but the data are insufficient to determine quantitative trends. 
displacement, 5, =(p- 
in moving from maximum to ‘residual stress—from Point A to Point B in Fig. 
3—is different than that for another pile diameter. The pile segment and adjoining 
- soil up to the point of soil failure, Point A, will be governed by ‘‘elastic’’ 
- conditions. Once soil failure is reached, the #-z curve of the pile segment will | 
be governed by movements along the failure surface and by “‘elastic’’ conditions 
for the soil away from the failure surfa ce. ‘This implies ~ that once failure | is 
reached and if the load with continued pile displacement, 
_ the soil beyond the failure surface will ‘‘snap back’’ a small amount as its 
_ elastic energy is released. As the pile diameter increases, this action increases 
_ the relative displacement along the failure surface and causes the soil for a_ 
large diameter pile to reach a residual level at a displacement less than that 
for a small diameter pile in the same & 
Pile displacement at the beginning point of ultimate load transfer, Point A, 
in a given soil increases with an increase in pile diameter, and the value of © 
: * is nearly constant or decreases, due to “‘snap back’’, with an increase in 
| - pile diameter (22). Therefore, » may fen with an increase in pile diameter. © 
These effects combine to result in pile length, or pile-soil stiffness, and | have 
greater impact on the capacity of larger diameter piles. t lo 
__ Typical values of » that we have observed from load tests range from about 
2-5. This range of » values can be combined with the range in & values to 
determine the impact of the softening phenomenon on axial capacity. This effect — 
is shown in Fig. 6 in terms of mobilized friction ratio and pile-soil stiffness 
- for an € of 0.75 and three values of p. _ These results show that the relative | 
"magnitude of reduction in the friction ratio decreases with an increase in pile-soil — 
_ stiffness ratio. For 7, greater than about 10, the additional friction ratio decrease a 
is less than 15% if p is less than about five. The above results show that Pet | 


*s and the shape of the f-z curve can explain at least part of the nag ae. 
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Lateral Pile Movements. movements installation can al also 
contribute to a length effect by affecting the magnitude Of fina max along the p pile 
length. The lateral vibrations result in low lateral stresses acting on the pile. 
‘The more pile that is driven, the greater the potential deterioration of the stresses, 
especially along the upper portion of the pile. The amount of lateral movement, 
the depth to which it extends, and the consequences are affected by the method | 


4 
_ _and details of installation—driven, jacked, alignment control, free standing length, : 


me ‘pilot hole, bending stiffness of the pile, soil strength, and soil stiffness. Any 
th effect due to lateral pile movements is not readily amenable to analysis, a 
but the effective stress approach to pile capacity that uses cavity expansion — 
to simulate pile installation has potential to make advances in this area. sapean, 
Comment.—The foregoing concepts can be used to explain the of 


q ‘MOBILIZED FRICTION 
025 0.50 0.75 vit 1.0 


PILE-SOIL STIFFNESS 


| 
in Theoretical Length Effect 


length on axial capacity, but unfortunately, data are unavailable to test = j 


hypotheses. Thus, we must resort to empirical n methods tempered with judgment 
more detailed soils data combined good load tests on 5 


Base. ‘to come compute the axial capacity of. 


4 


correlated to both pile and pile- stiffness, 7,. Linear regression 
analyses were used to determine the most effective combination of dependent — 
a, B, or A and independent parameters L, , In (L), Ts; or In 7, for pile capacity — 
_ predictions. . The piles analyzed in this study were loaded to failure within a i 
few hours after initiating the load test. They had diameters between about 6 
in. (152 mm) and 30 in. (762 mm), lengths between 8 ft (2.4 m) and 333 ft 
%. (102 m), and included pipe piles, concrete piles, and timber piles. The available _ 
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IG. 7.—a Profile (1 ft = 0.305 
load test iain in include pipe pile data from Vijayvergiya and Focht (32), timber 
e and concrete pile data from Flaate and Selnes (15), and data from Meyerhof . q 
and Murdock (25), Ito and Koizumi (18), Hutchinson and Jensen (17), Stermac, 
and Devata Tomlinson Danys Muktanbhant, and 


sumber of good load test data is relatively small compared with 
the number needed to in a significant manner the effect 


_ number of additional data points will not statistically af fect the regression orn 


> The period between installation and testing for the pipe pile data was’ “- 
excess of 30 days with the exception of the data from 44 of the 105 tests. 
It was less than 25 days for 30 of these 44 tests. Once a pile is tested to - 
_ failure, its load carrying ability may decrease to a residual capacity with an 
_ increase in pile movement. Subsequent testing to failure in this case can result . 
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FIG. Profile (1ft=0.305m) 
in peak loads less ‘ieee initial peak and closer to the residual capacity, 
as evidenced by data from Meyerhof and Murdock (25), Bergemann (3), and = 
ecm _ Ifa pile is redriven to a deeper penetration after long delays, the peak capacity 
; a may be less than could be developed if the pile were originally driven to _ 
_ deeper penetration with minimal delays. Such behavior tends to occur when _ 
the soil exhibits a strain softening behavior and the pile is subjected to a 


movements. Therefore, piles loaded tof failure and Seen driven to dongs 
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penetrations for later ony: capacities the deeper 
$ penetrations th than would be developed if t the piles were ‘originally driven to ; 2 
- the deeper penetration. If load test data are excluded for f piles with less than 
= set-up, piles having oversized plates at the tip, and piles with need a 

-_testing-redriving-retesting, the data base would not be sufficient to lend statistical | Pa 
"significance to any empirical procedure. Therefore, the influences of set-up 
time, oversized closure plate, and loading history had to be neglected ~~ 

regression analyses, in 1 recognition of fact that the results will 
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FIG. 9.—) Profile (1 ft = 0.305 

Regression Analyses. _—The load test data are Presented in Figs. 7-93 in terms 

oy of depth profiles of average adhesion ratio (a = ratio of average adhesion, 

f,,, to average undrained shear strength, s,,); B values = f/G,,,; \ values, 
and in Fig. 10 as a profile of \ versus 7,, the pile-soil stiffness term. The 
amount of set-up time, piles with oversized closure plates, and piles redriven 
to deeper penetrations or reloaded, are identified for the data shown in Figs. 
7-10. To compute ™;, we to took Snax 28 the mean undrained shear strength and 


as 0.1 in. (2.54 as sufficient data were not available to determine 
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. The value ot of 0.1 in. (2.54 mm) is aameniibia the pile — ; 
_ at maximum load transfer used with empirical t-z curves (8, 9). The elastic moduli 
_ of the piles were taken as 29,000,000 psi (200,000,000 kN / m7”) for steel, 2,600,000 3 
(18,000,000 kN/m/7) for concrete, and 1,760,000 psi (12,000,000 kKN/m7) 

d for timber. The wall thickness of some pipe piles w was not known, thus 7 > | 
not be ‘computed for all tests. nil Yo af ative od! " 
The data were separated according to pile type and soil stress history. No 
statistical difference in pile capacity factors (a, B, 4) was found between the — 
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different pile types for the data in this although others have 
found timber piles to wae higher unit friction values than steel piles (4). 
A value of the average s,,,,/¢,,, of 0.4 was arbitrarily taken as the dividing = 


were print to a, to pile ‘penetration, 
in of pile penetration, pile-soil stiffness, 35 and In of 7. In each case, two 


greater than 0. 4. In much of this data, the 


7 
04 and one for s,,,, 


- The rate of decrease in the pile capacity coefficient (a, B, A) is expected 
to decrease with | an increase in the independent parameter (Fig. 3). _ Therefore, * 


even ‘if the regression results in terms of 7 7, OF “penetration are e slightly more 

Statistically significant than those in terms of the logarithm of these independent 

_ parameters. The arithmetric values of the independent parameters only provided 

a better correlation for the normally consolidated soils and then the differences 

were small. . Therefore, the results of t the -Tegression analyses are summarized — 
= 2 for the logarithms of the independent parameters. Table 2 


lation 


coeffi- 


= 0,296-0.032 In (L)" 
B = 0.468-0.052 In(L)” 
= 1.486-0.126 In (L) 
= 0.178-0.016 In 7, 
= 0.278-0.028 In 
a = 1.024-0.070 In 7, 
d = 0.488-0.078 In (L) 
a = 1.012-0.103 In (L) 
= 0.232-0.032 In 7, 
= 0.685-0.061 In, 


“Pile L,infeet. 

maturallogarithm. 
also the correlation coefficients of the regression analy 


‘The correlation coefficient provides a measure of the degree of fit for the 


_ assumed correlation function. Low absolute values of the correlation coefficient _ 


; _ the data or, if the standard deviation about the regression line (scatter) is small, 

that the dependent parameter is not significantly affected by the independent 

_ parameter. Standard deviations of the different methods cannot be directly i 
compared as the magnitudes of the means of a, and differ by almost 
a magnitude. Therefore, the ratios of predicted to measured values a, B, and 
A were computed. The mean and standard deviations of these ratios are also 
in Table Both the deviation of the ratio of 
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| be 
j 4 data Mean | tion 
| 1.06 | 0.23 | 1.29 
058 | 1.05) | 022 | 1.28 
nF | 110 034 | 1.44 
4 4 47 70.68 | 1.05 | 0.24 
—-065 | 105 | 030 | 
4 | 105 | 0.25 
48 1.09 | 034 _ 
— on 
- = be interpreted to mean that either the regression function does not represent _ i 7 
if 


computed t to o measured ‘capacity must be be. considered iia evaluating | the confi- 
: S dence level of a predictive procedure. When the mean exceeds one, the sum — 
_ of the mean and standard deviation provides a measure by which to mcr 
the different methots: the The follow the s sum, the better the method. 


«FAG. 11. —Comparison of Regression Curves (1 ft = 0. 305 m) | 


‘emtalee that are slightly better than the a method as reflected by the correlation a 
_ coefficients and the sum of the mean and standard deviations of | the ratio — 


3. For overconsoiidated soils, the a and methods provide 


C b. 4. The use of a soil- -pile stiffness ratio as an independent parameter provides > 


a ‘slightly better correlation than pile penetration for normally consolidated soils. 


B values are highly” sensitive to OCR, and this sensitivity contributes 
to a low correlation between B and either pile penetration or 7,. The A value 
_ also influenced by OCR, but the influence of OCR is less, as as the a method 
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magnitudes of a, and are not highly sensitive to soil strength for 
normally consolidated soils. Therefore, results from a correlation of these 


parameters with or 7,, in consolidated soils, should 


™ uses both undrained shear strength and initial effective stress in the correlation; =§ [7m 

=the and a methods use only one of these two 

TABLE 3.—Comparison of Methods 

Diam- | thick- | pounds | tra; | 

| 4& | 10 | 300 | 79.3 0.108 | 0.96 
(84 a $3.7 | 0.100" | 0.114 

| | | 284 | 01200 | oie | 10 

- —— 


_ between measured £ values and only pile penetration or 7, in overcoisolidated 
soils was found to be low, and the results are not included here. The load 
- test data for overconsolidated clays are for short piles, and the | targe variations | 


Less data scatter and greater ¢ a values are expected for test results from long 


piles in deep, overconsolidated profiles. §|. | 


The regression curves for \ in terms of pile penetration and pile- soil stiffness _ 
are compared in Fig. 11 and 12. The regression curves provide higher \ values" a 
‘ for overconsolidated conditions than normally consolidated conditions for shallow | 
- penetrations, but at a pile penetration of 70 ft (21 m) or a 7. of 24, the regression ‘| 
curves for: normally consolidated and overconsolidated cross. This crossing 
is due in part to differences in the distribution of data points between the Jf 
two conditions for s,,, /G,,, less or greater than 0.4. 
_ The influence of the A correlations with pile penetration and 7, on computed 
pile capacity was examined for five offshore conditions | encompassing the 
nig og encountered range. ‘These « examples, which are summarized in Table _ 
3, show that 7, for offshore piles ranges from about 5-80. A value of ur 
& (2.54 mm) for u* was used to compute 7,. The ratio of average shaft 
friction determined from the 7, correlation (Case 4) to that in terms of pile 
_ penetration (Case 1) ranges from about 0.95-1.15. From statistical ee. 


(t-z support the observation that a length 
an exists on the friction capacity of piles in clay. The magnitude of the 
effect influence is affected by soil stress-strain behavior (€, pile-soil 
‘stiffness | (m3), lateral pile movements during installation, overconsolidation ratio, 
and other factors thought to be secondary. Parameter analyses demonstrate 
that average friction for a pile is not likely to be related to a single parameter 
representing pile-soil stiffness and soil stress-strain characteristics. Unfortunately, — 
- sufficient data are not available to test a theoretical model. Therefore, empirical - 
approaches based on theoretical concepts must be used to evaluate pile capacity 
Re A statistical analysis of the available data show that the \ concept, either 
in terms of pile penetration or 7,, provides the most consistent and reliable 
single method for computing the axial capacity of piles in both normally 5 
consolidated and overconsolidated soils. The computed capacity is for loading 
‘rates where the failure load is reached within a few hours after initiating the © 
a Until additional pile load test data become available in sufficient: quantity 
- to support alternative methods, we recommend that the axial static capacity 
of long” piles driven into normally consolidated and lightly aaa | 
(Sim /F vm < 0. 4) clay formations be computed with the Lambda or B correlation — 
with In T; ‘(Cases » 4 or 5, Table 2) that provides the greater er capacity. For piles a 
in heavily overconsolidated clays, we recommend using either the Lambda or 
a correlation with In Ts (Cases 9 or 10, ‘Table 2). The value of 7, 


[ . 7 a provide an indication of whether or not there is a length effect. The magnitude 
of B is influenced by the overconsolidation rati able herefore aco 


GT11 
computed in the same way as was dom. for the correlation study. ~<a 
ecause of short set-up times, oversized closure plates, reloading, and redriving 
some of the | test. piles, "especially those with large penetrations and large 
ues of 7,, ,, the regression curves tend to “underestimate the actual capacity —§ 
long piles without closure plates and with long set-up. Therefore, we further - 
recommend, based on personal judgment and the available data, that the minimum a 
‘A value be limited to 0.14 and the minimum f value be limited to 0.23. " 
_ Fig. 13 shows a comparison of B profiles in normally consolidated clays for 
e 1, with limited to 0.14, Case 2, B limited to 0 0.23 1.23, and the Flaate < 


‘ot 


= 03 AND 
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and Selnes method. A value of 0.3 for B, (s, fe. ,) was used to generate the 
B curve for the Case | curves. Three methods (Cases | and 2 with limiting 
values and Flaate and Selnes) are in reasonably good agreement. 
_ The recommendations in the preceding paragraph apply only to piles whines 
"oversized closure plates and piles with full set-up. The minimum value of s,,, /¢,,,, 
from the load test data was about 0.2, thus the recommended procedure may > 
not be applicable to piles in soils with lower ratios of s,,, y . The influence — 
rate and loading alterations to the capacities 
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and with « loading. due to degradation of the shear ‘resistance 


_ The siutiiin data base of good quality pile load tests is very limited and 


- design of marine structures. Typical of the 1 range in pile and soil conditions 
are pile lengths of 200 ft-600 ft (61 m—183 m), pile diameters of 3 ft (0.9 m)-12 | 

- ft (3.7 m), pile wall thicknesses of 0.5—3 in. (13-76 mm), soils in underconsolidated 
i to heavily overconsolidated states, and soils with shear strengths of 0.05- Is 
ksf (2.4-718 kN/m’). The cyclic nature and rate of wave loading experienced 

is by piles supporting offshore Structures are a stro1 strong contrast to the 1 monotonic, a 
slo slow ' loading rate on most test piles. Thus, the ge geotechnical engineer has been 
and will continue to be faced with the challenge of computing axial pile capacities 


for conditions for which there is little or no prior load test experience. 


computed and 


* ‘Soil variability and disturbance of the soil surrounding the pile | 

Distribution and magnitude of ultimate skin friction values 

Changes in the soil properties with time (consolidation, M: 

thixotropic and seasonal effects, and load history) My 
D Distribution and magnitude of residual stresses in the pie t 

Pile installation details (displacement characteristics, time delays, 


___Sequence of installation) 


“Actual capacity may differ by +10% from computed capacity, 
ii Some of the major factors that contribute to differences between actual and 
aa capacities of static monotonically loaded piles driven into clay ¥ 


in evaluating pile load test data are summarized in Table 4. Subjective estimates, 
~ which reflect our opinion of the present state-of-the-art, . are shown in Table 
4 for the relative uncertainty of each of these factors on pile capacity for 


and soil conditions where we have previous experience and where the 

piles are installed with good construction practices. Considering these factors _ 

and limitations of the analytic models, we expect measured Capacities in most 


cases to be between 0.7 and 1.3 times predicted capacities. 
_ Some methods may be more appropriate for certain soil conditions than pn 


The difficulty lay in establishing which method is the most appropriate for — 
>. which conditions. The first step is. to recognize the limitations of the available 
‘pile capacity methods. If a method is not appropriate or the uncertainty is 
7 unacceptably large, load tests may be needed to guide and confirm design. 
wan load tests are not practical, the geotechnical engineer must make rational _ 


| | 
40.1 
| 
| 


decisions | about pile | most ost appropriate method | 


on his experience and analytic studies. Although recent research on effective | 

, stress methods with cavity expansion and reconsolidation offer promise, they _ 
are not yet ready for use as a standard design tool. Nevertheless, these methods _ 
can be a useful analytic tool to gain insight to aide the decision- -making process. — 


ures and increase the reliability and ‘cost: - effectiveness of marine pile 
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is an expedient means of computing the axial movement of a i under pow 
load. Procedures are available to generate the relationships between shear stress 
at the pile shaft (load transfer, 1) and pile displacement, Zz along the pile shaft — 
and at the tip (3,4,6,8,9,10,13,16,21,22). The most commonly used procedures, 
_ however, are empirical and based on data from tests on short land piles, usually 
_ less than 100 ft (30 m) long, with diameters less than 18 in. (0.5 m). Pile diameter, 
- axial pile stiffness, pile length, and distribution of soil strength and stiffness 
along the pile are all factors that influence the t-z behavior. Data from instrumented — 
pile tests do not span a wide enough variety of conditions | to establish au 
comprebensive data base for all the variables in the theory. _ ae Rn) 
_ The theoretical concepts presented in this paper were developed to provide 
a basis for f-z criterion that would be applicable to a variety of pile and soil ; 
conditions. The validity of the tested with case 


DEVELOPMENT oF f-z Curves" 
_ General Comments.—The at displacement of point 
= are influenced by the stress at every other point; the degree of influence 
7 = decreases with distance from the point of interest. The subgrade-reaction approach 
assumes that the displacement at any point depends only on the stress at that 
point. Soil-pile interaction is represented by springs that relate shear’ stress 
on the soil-pile interface and displacement of the pile (Fig. 1), 
_ Success in developing realistic t-z relationships for a pile depends on the 
accuracy of the ultimate load-transfer values of the soil (pile capacity), the 
Getcibution of those values along the pile, , and the displacement characteristics 
'Mgr., Special Projects Group, McClelland Engrs., Inc., Houston, Tex. 77081. 
?Grad. Research Asst., Univ. of Michigan, Ann Arbor, Mich., formerly Geotechnical 
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of the American Society of Civil Engineers, ©ASCE, Vol. 1 107, No. GT11, November, 
1981. ISSN 0093-6405 /81 /0011-1543 /$01.00. 
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of the soil during load transfer. For this ‘study, ‘oii techniques were used 

to predict the ultimate values; therefore, only the load- transfer characteristics 
| 

as they relate t to are re discussed i in n detail. 


cal model on elasticity. Postfailure relationships were modeled considering 
the residual stress- deformation behavior at the soil-pile interface. i 
 Prefailure t-z Curves.—The theoretical formulation discussed below uses a 
concentric cylinder approach with the following key re and ame 


as described Randolph and Wroth (18): 
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FIG. 1.—Finite Difference Model of Axially Load ed Pile 


. The displacement pattern of soil can be modeled as concentric ae 


2. Radial soil displacements due to pile loads are assumed negligible when 


> to vertical soil deformations. Therefore, simple shear condition prevails 

om stress at distance r; +, = the shear stress on the piles -soil idiathiaee 


Shear stresses are negligible beyond radial distance (or zone of 


_ influence), and the soil does not deform beyond | that point. ~~ ati 


| 
| 
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‘on the pile; in which / = pile length; v = Poisson’s ratio of soil; and p = i 


atio of the soil shear moduli at depths //2 andthe pile tip. = 


in w which z, of of shaft element. G be a function 
‘ distance as a result of disturbances caused by pile installation and Sa 
soil behavior. Both of these possibilities will be discussed in later sections. 


soe ag failure, the slip planes that develop invalidate Eq. I. Up to the point 
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_ 
FIG. 2. tate Modeled as Shearing g of C 


Be = 


te) 


Eg. 2 shows that the displacement is not only a a function st a 
- soil stiffness, but also pile radius and length (indirectly by r,,). These eset 
not directly accounted for in empirical f-z criteria. poe 
i wa Variation of Shear Modulus with Depth. —tThe soil moduhes varies along the 
wt pile, with radial distance from the pile, and with pile load as a result of natural ‘ 
a variations, soil disturbance during pile installation, soil consolidation during ~ 
setup, and nonlinear stress-strain response. The f-z response at a particular 
depth, however, is taken to be controlled by soil modulus conditions at = 
depth. Using this assumption, Randolph and Wroth (18) approximated a linear 
_ increasing variation in soil modulus along the pile shaft. The load- displacement _ 
_ response computed was in good agreement with more rigorous solutions. 


Therefore, this basic in the in approach should account for vertical 


Tain can De modeled Eq. 1. 


* Radial Variation of Shear Modulus. —A . typical distribution of the soil modulus | 


is shown in Fig. 3 for four conditions: (1) Immediately after installation for 
=. loads that result in small soil strains; Q) immediately after installation for _ 


— 


AFTER 


RADIAL DISTANCE 
RADIUS 


— SMALL Loans 


i 


MODULUS RATIO, 


—Variation i in Soil Modulus Ratio 

‘pile loads near failure that result in large soil ae (3) ) after c saaidiiiiainn _ 
for small pile loads; and (4) after consolidation for pile loads near failure. = 
‘The variation in modulus due to to pile installation and | soil: seemed can 


_ 
| 
| 
8 if 
— 
me 
= NEAR FAILURE LOAD 
— | 
~ NEAR FAILURE LOAD 
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consolidation show that the size of the zone of (large stress changes) 
due to installation extends to about 20 pile radii from the pile, but soil failure 
extends only to about 4 pile radii. After consolidation, the shear strength next 
to the pile, which is predicted with the consolidation analyses, may exceed © 
_ the original value with the strength decreasing away from the pile cs ne 
4 If the soil modulus is proportional to the undrained shear strength, then the 


distribution of the ratio of the shear ‘strength after installation to initial shear 
seen can be used to adjust the initial in situ modulus at small strains to e 


om 


several 


ised wit 
ath Eq: DISTANCE FROM 
PILE RADIUS 


|.—Idealized Radial Distribution of Soil Maduius Ratio 


an average modulus ‘that will provide the same soil displacement, z,, next to 
the pile wall as a solution for the actual distribution. For a linear variation 2 
i 


in modulus (Fig. 4) the average soil modulus i dulusis 


duly 


im in ‘which GC... = ‘equivalent shear 1 modulus for the pile-soil system; G, = shear _ 
modulus a great distance from pile; and r, and M, are defined in Fig. ry The | 
Pie lily variation i in 1 modulus i is not known witha high Gegres of of precision. Therefore, 


+1 ) 


— 
g _ be estimated from results of cavity expansion to simulate changes in soil stresses : a 
= during installation and consolidation theory to simulate the stress changes due a 
7 = to dissipation of pore pressures induced during installation. Available solutions Jf 
q 
at 
time after installation. T 
oF 
knowledge of the stress-strain behavior of the soil. 
| baw how gy 


linear approximation of this variation is for most cases. 
As an example, equivalent, _ weighted moduli were determined with Eq. 3 
for a stiff, overconsolidated soil with a sensitivity y of 1.5 and a maximum gain 
in strength after consolidation of 30%. The approximations to the radial variations a 
of modulus ratio are shown in Fig. 5 for conditions after consolidation. These 
estimates were made using cavity expansion and consolidation results by 
Randolph, gprs and Wroth (17) and a nonlinear stress-strain curve. Values — 
‘M,, and Gye / Gu at small and large strains are shown in Table 

TABLE 1. —Radial Variation of Soil Modulus Ratio Parameters® and Average Values 


_ Immediately after installation 


Large load, f = 0.67s, 


Small load 


AFTER LOADING 04 sy 
AFTER LOADING = 10 sy ' 


XIMAT 


PILE RADIUS 


FIG. 5. —Estimated Radial Distribution ot Soil Modu 


For conditions ad consolidation, the nonlinear effects are more important 
than the installation effects. Therefore, further assessment and consideration 
of the effect due to installation ' will be neglected in the following work, and 


7 _ Nonlinear Stress-Strain. —The stress-strain behavior of some soils can be 


of 


max 


| | | 
| | 
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in which G: = secant shear ‘modulus at an stress ; G, = initial 

_ shear modulus at small Strains; R,= stress-strain curve-fitting Constant; and 

a = shear stress at failure. By combining Eqs. 4, |, 4, 1, and the distribution of of — 


with |r, the equation for t-z response becomes 

g which o> =7,R, Tenax: A similar "equation: be for a 


Eq. 5 can be used t to generate prefailure t-z curves. “Although 
neglects radial heterogeneity caused by the installation process, these effects, a 
7 after a long setup for driven piles, may be minor compared to en we 


in the properties of the undisturbed soil and to the approximations used in 


and stress-strain nonlinearity share equal importance, both effects can be 

combined as discussed previously. 1 

Post-Failure Curves.—We must t address several related problems to develop 

an analytic model that includes t-z behavior at and beyond failure. These problems 

_ include: (1) Maximum shaft friction; (2) pile displacement (or strain) at which — 

the maximum friction is mobilized; (3) residual shaft friction at large pile 
displacements; and (4) displacement behavior between maximum and residual 7 

stresses. In this study, we computed the maximum shaft friction by conventional 

means and the pile displacement, z,,,,, at which the maximum friction is mobil- 


Some conventional ‘methods to compute shaft friction | provide 
. These methods account empirically for 
differences between maximum and residual load transfer and the degree of 
shear resistance between these values that is mobilized at maximum load on — 

‘J the pile. Also, these conventional methods provide average resistances along 

4 the shaft that may not reflect the actual distribution along the shaft. The 

distribution of shaft friction along the pile is important to the load- “displacement — 

- behavior of the pile. For example, a shaft friction of 1 kip/sq ft (48 kN/m? ) 
may be used to compute the capacity of a pile in a stiff Beaumont clay. The 
distribution along the shaft, however, may be zero at the ground surface and — 
increase linearly to 2 kip/sq ft (96 kN/m’) near the pile tip. Some discretion 
should be used to establish the distribution of shaft friction along the pile ' when 

4 using conventional methods to compute the maximum shaft friction. 
Once the soil reaches a state of incipient failure, additional deformation at _ 
the pile wall results from a shear strain in the soil with a rotation of principal 
planes, movements along slip planes, or a combination of these conditions. 

_ When slip planes develop, the soil behavior cannot be defined in terms of 

, ‘stress and strain. Therefore, Eq. 1 cannot be used to define the t-z response 

_ for postfailure conditions when slip planes develop. One approach to postfailure 

7 ~ behavior is to model a section of the pile-soil system in a direct shear test _ 

-s some other laboratory simulation. This is done by replacing one half of - 

- the soil sample normally used for the test with a representative piece of steel, 
—_— or wood. ‘Ifthe failure is in the soil rather than at the soil- _ interface, 
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a ‘conventional direct shear test may to define the load 
transfer- displacment response for post-failure ‘conditions. 
. Direct Shear Tests on Sand and Concrete.— Desai and Holloway (5) and Kulhawy 
and Peterson (11) present load-displacement curves for sand and sand-concrete 
interfaces. The displacements at which the residual shear stress was reached © 
are typically between 0.035 in. (0.9 mm) and 0.05 in. (1.3 mm) beyond the 
haan at maximum shear stress. While the data from these studies « are 


by no means inclusive, the data indicate that residual shear resistance i ai sand 
range between about 0.80-1.0 times the maximum shear strength. = 

Direct Shear Tests—Clay.—Residual resistances for undisturbed range” 

from less than 60%-100% of the maximum resistance (12). The magnitude of 

_ strain softening of remolded material may differ from that of undisturbed material. 


—— 


“The degree of softening may also be influenced by test type (simple s shear, 


“ELASTIC” 


outw od) 6.—Exaggerated isp ro ie od 

As part of this study, several consolidated-drained direct shear tests were 

performed on undisturbed and remolded samples of highly plastic clay. We 

tested a clay-clay interface instead of a clay-pile interface because several model 
and full-scale load tests have shown that the failure surface for piles in clay 
- occurs within the clay, a small distance from the clay-pile interface about 20% q 

_ Results of the direct-shear tests indicate a strength decrease of about 20% 

Pe peak to residual for this clay. The deformation necessary to develop the 

_ residual value after reaching the maximum stress was about 0.1 in. (2.5 mm). — : 
- Interpretation of Direct Shear Tests.—The direct shear test is not a perfect 

"simulation of the load transfer- displacement behavior of a pile segment. The 
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boxes also result in differences between the simulation and prototype. Further- 

more, the stress-displacement curve generated from a direct shear test simulates — 

the condition along a failure plane (Fig. 6), and does not account for elastic — 

strains that occur in the soil some distance away from the pile. eens bas : ss 

_ The displacement, 8,, that occurs in moving from the maximum to ‘residual 
stress in a direct shear test can be compared to the f-z displacement, Az, that 

occurs in moving from oe maximum to residual load- cones wees for a 


- surface, and by Eq. 5 for soil away from the failure surface. This implies 

that, once failure is reached and if the load transfer decreases with continued — 

SCALE FROM DIRECT 
SHEAR TEST RESULTS 


DIRECT SHEAR TEST 


SHEAR 
STRESS, T/T 


‘PILE 
_ pile displacement, the soil beyond the failure surface will ‘ “snap back” a small 
amount as its elastic energy is released. This action inevenste the relative 
_ displacement along the failure surface and causes the soil to reach a residual 


level at a a pile less" than ‘that | given by « direct results. The 


- which Az, = the elastic rebound, which is cunpuned from Eq. 5 , taking the 

“¥ The maximum difference between Az and 8, is less than (10% of 8, for a 
12-in. (0.3-m) diameter pile and typical soil conditions. The value of Az, increases 
with an increase in pile diameter, and 8, is probably insensitive to pile diameter. it 
Therefore, for pile diameters greater than about 18 in. (0.5 m), the magnitude _ 

of Az, should be evaluated to determine if a correction is needed. Otherwise, - 


THEORETICAL 
— total normal stress is held constant during the direct shear test. This stress fa i 
2 boundary condition is different from the displacement boundary condition (nearly * 
«Zero radial movement) during pile loading. The stress concentrations caused nd 7 
ne ma Onta area he lahorato and nacing he Pen hear 
q 
q 
; a the point of failure, prescribed by Eq. 5. Once failure is reached, the f-z curve ' 
he ni seoment is soverned by th irect sh ear tect results the fai 
= value OF 0, Irom direct shear tests 
| 
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| Theoretical te 

- in two parts: prefailure and postfailure. The basic concept is summarized in 

Fig. 7. Prefailure #-z is governed by Eq. 5. This equation is flexible and | 

= be altered to simulate changes of ‘Shear modulus with shear ‘Strain and z 


“used for square or H- -piles when us using 1g Eq. 5. The load- -displacement aintiahins 
_ require a maximum unit load transfer value, f,,,,; shear strength, s, ; soil modulus, — 
Gag and stress-strain constant, R,. The value of ¢,,,, may be less than 
for example, for piles in stiff, overconsolidated 
f tmax 1S COMputed in the same way as it would be for shaft friction in pile oo 
capacity computations, but some discretion should be used to select the distribu- — 
tion of ¢,,,., along the pile. For example, a a linear increasing distribution of aa 7» b 
is probably more appropriate for short, driven piles in stiff, overconsolidated _ 
clays than a uniform value, even though the shear strength may be constant ~ 
over the pile length. The average of the distribution should be consistent with — - 
the value obtained from conventional pile ‘capacity procedures. If laboratory on 
data show that the residual shear resistance is less than bax» the pile capacity 
ee is used to provide a value of shear resistance between a, Ya ‘ 
Thus, the value of ¢,,,, should exceed the value obtained from the pile capacity : 
procedures. Until more refined guidelines and procedures are developed, oe “a 
value from the pile capacity procedure is used to represent the mean of i 
The difference between the capacity generated from the t- -z analysis, — 


res* 


_ Trial and error can be used to revise the estimates, if necessary, until the 
for the two methods are within 


rigid o1 or very compressible piles and for large between and 


(7). data c can be used to determine the: nonlinear response, as reflected 
"Ones the failure stress has been reached, the post-failure behavior can be 
“estimated from results of direct shear tests. Limited data incicate that the 
deformation that occurs in moving from “ie to ha is between 0.03 in. ©. 8 
_ mm) and 0.05 in. (1.3 mm) for sand and about 0.1 in. (2.5 5 mm) for clay. ‘ee 
Comparison between Theoretical and Empirical t-z Curves.—The deformation 
needed to develop full shaft friction, z,,,., varies considerably. This variation 
can be accounted for by the soil shear modulus, G, and the stress-strain param- 
eter, R,, pile length, and pile diameter, used in the theoretical approach. None 
of these factors are directly considered in empirical formulations. » 3 ee 
‘Three examples of t-z construction (one each for sand; stiff clay, and soft — 
» are shown in Fig. 8. For these examples, ‘he Post- peak response was ~ 
approximated by a straight line from ¢,,,,-90% of the maximum for pile movements y 
a than z,,,. The curves for a sand show that empirical procedures i 4 


result in a much softer f-z spring than indicated by the theory presented here 
on which the empirical procedures were 6-in. -mm)- 18-in. (457- 


The curves for clay are within the range of the aga ay curves, but the | theoretical 
mm) Coma — with less than 100- ft (30-m) penetration. The comparisons : 
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"THEORETICAL CURVES © 


z Response at Pile Tip.— —The Q-z curves at the tip used j in the ‘subgrade 
= approach have been developed from instrumented pile tests and labora- _ 
e tory model tests. While the results of such tests yield useful information on _ 
7 ae load-settlement performance, a theoretical approach is needed to study | 

the fundamental factors that influence the behavior of Q-z response and to_ 
_ apply the Q-z concepts to conditions beyond the data base on which the the empirical — 
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«AIG. 8. ammo of t-z Curves (1 in. = 25.4 mm; 1 ksf = 47.9 kN/m apt + 


‘Thee elastic | solution for a a rigid punch (20) car can 1 be ‘used to model Q-z oa 

tip diameter; E = = Young’s modulus of 


— 
at VIJAYVERGIVA | — 
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L = influence coefficient. _ An appropriate secant modulus 


is used in Eq. 7 to determine a value of z. For a rigid punch on an elastic 
7 half-space, the value | of J, = 0.78. Studies reported by Randolph and Wroth 
(18), in addition to those by Vesic (21), show that J, is reduced, due to depth | 
effects, to values typically between 0.5 and 0.78. ee were: Ae 


_ To verify the theoretical t-z approach | and the procedures to select soil input i 
parameters, we predicted load-displacement response for four single piles tested 4 
_ in a medium dense sand at a site on the Gulf coast, and single piles in stiff — 2 


at the of Houston Hendon, » The soil and 


FIG. -9.—Soil and Pile” Conditions at Houston Test Site (1 ft = 0. 305 m; 1 kip/oq 


of predicted and measured load-displacement curves for each test. All the case 
_ studies were for short-term loading. The computer program AXPILE was used © 


Sand 

4 Pile and Soll Conditions.— —Load tests on four Armco Hel-Cor pil Cor ibis were. 
performed in 1969 at a Houston site. Three of the 12.25-in. (311-mm) pean 
piles (Piles 2, 4, and 5) were driven to 57.5-ft (17.5-m) penetration, and a fourth | 
pile (Pile me was driven to 46- ft (+m) penctration. The piles had an Sone 


TOF cach OF these tests are described briefly along with a 
{ 


: 
THEORETICAL C CURVES. 


‘The eolta at the site consist primarily of a loose- to-medinndease silty sand. rf 


Profiles of standard penetration blow counts are shown in Fig. had for the two a 


‘with s an » earth pressure coefficient of 0.7 and a friction angle of 30°, except 5 
_ for depths greater than 41 ft (12.5 m) for Piles 2, 4, and 5 and for depths _ 
"greater than 34.5 ft (10.5 m) for Pile 3 for which a friction angle of 35° was 
used. An end bearing coefficient of 40 was used to compute the tip resistance, J 
which was about 50% of the computed capacity. 
_ The f-z curves were generated with Eq. 5 using an R, of 0.9 and the soil 
modulus profile (Fig. 9). Any loss of shear resistance for pile displacements _ 
21) 


TEST PILE 


MOVEMENT, IN. 


FIG. 10.—Load- -Displacement Piles at Site (1 in. = 25 4 


mm; 1 ki 4.45 kN me = 

greater than z,,, was arbitrarily neglected. The load-displacement (Q-z) curve © 
at the pile tip was computed using Eq. 7 with E = E,,,, for values of Q up 
to 0.25 E=05£,,, for = 0.25 E,,,, for and 
straight line segments at intermediate points on the Q-z curve. A Poisson’s 
ratio > of 0. 3 was used to convert G to bivaw it 


7 


curves are ¢ compared (Fig. 10). The ‘computed results compare well with the: 


| | 


in 1979 and ‘1980 at a site on as "University of Houston Campus to 
the performance of single piles and groups of piles driven into stiff, overconsoli- 
: dated clays (14). Eleven 10.75-in. (279-mm) diameter closed end pipe piles were 
driven to 43-ft (13.1-m) penetration; nine in a square group with a 2.75-ft (0.84-m) 
- center-to-center spacing; and two single piles located about 12 ft (3.7 m) from 
the center of the group and on opposite sides of the group. The geometry 


of the pile group and pile roperties are Shown in Fig. Us Details of the 
pus p 


repeated here. The piles were incrementally loaded at about 20, 82, and 110° 
days after installation. Measured and predicted results are compared for the 
ey soils at the test site consist of strong clays and sandy clays. The — 


_ overconsolidation ratio decreases with depth, but exceeds 3 even at 50-ft (15.3-m) 
depth. The interpreted modulus profile, which is shown in Fig. 11, was based a 


INTERPRETED UNDRAINED INTERPRETED ELASTIC 
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PIPE PILES 
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on 
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FIG. " .—Soil Conditions at University of Houston Test Site (1 ft = 0. 305 m: 1 kip/sq 


he 
on cross-hole reed The strength profile was based primarily 
on unconsolidated-undrained tests, as the method used to predict pile capacity 
is based on this kind of strength test. The cone data, however, revealed a 
trend in the strength profile different from the laboratory data, as shown in 
_ Fig. 11. This difference would not affect the estimate of shaft resistance with | 
the procedure used to compute pile capacity, but the difference would affect 
_ Anundrained shear strength of 6 kips/sq ft (287 kN / m7’) was used to compute 
- the ultimate end bearing resistance. The average shaft friction was taken from 
local experience as | kips/sq ft (48 kN/m7’). The shaft friction was assumed 


_ to increase linearly with depth from zero at the ground surface. = © _ 


| 
| 
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FAB. 12.—Load- Responses for Single Piles at of Houston 
To in. = 25.4 mm; 1 kip = 4.45 kN) 
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‘The t-z curves were generated with Eq. 5 using an R, of 09 9 and “the soil 


modulus profile in Fig. . 11. The development of the t-z carves arbitrarily neglected — 
any loss of shear resistance for pile displacements greater than z,,,.- 

load-displacement (Q-z) curve at the pile tip was taken as a binlinear curve ~ 
using Eq. 7 with = and E = E,,,, to compute Poisson’s ratio 

Load-Displacement Curves.—Fig. 12 compares the measured load-displacement 
_ responses for the two single piles with the computed response. The measured 

and predicted responses at working loads are in excellent agreement. bien paatens 

ultimate load is between the measured ultimate 


- Pile and Soil Conditions.—Load tests were made by Cooke, Price, and Tarr 
— Q) on 6.6-in. (168-mm) diameter piles jacked 15 ft (4.6 m) into | a stiff London 


clay at a site in Hendon, North London. Single “piles and a 1 pile group were 
tested. The ang in the group were in a line at a center-to- center spacing of 7 


& MEASURED. 
@ MEASURED. PLE 


HENDON 
sok 


MOVEMENT, IN 


ement for Single Pile 


to failure, but the Disctains hots at working loads. Details of the instrumentation, vi 
installation, and testing are given by Cooke, Price, 9: lillie 
The soil at the test site consists of stiff, overconsolidated London clay. 
“Interpreted profiles of undrained shear ‘strength and soil modulus at small strains 
are summarized in Fig. 13. The soil modulus profile was estimated with the 
empirical equations presented by Hardin and Drnevich (7). + ple cayeciny 
The ultimate pile capacity was computed using s, of 1.6 kips/sq ft (77 kN/m* ) 
at the pile tip, an average shaft friction of 0.94 kips/sq ft (45 kN/m7’), and 
_adistribution of shaft friction along the pile proportional to the undrained strength © , 
profile. The average shaft friction and its distribution were based on the load 
test data. The procedures used to compute the /-z curves followed those described 
for the University of Houston site, 
_ Lead-Displacement Curves.—Fig. 14 compares the load-displacement response 
for three piles tested individually to working loads. The predicted and measured 


results were in excellent agreement for working load conditions. The measured 


L 
{ 
ll 
FIG. 14.—Load-Displac Site (1 in. = 25.4 
| 
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the loading history and installation details between the two piles. Overall, the a 
predicted and measured results were in good agreement. 
Pile Group Tests.— —Analyses were also made to compare computed and 
measured response of pile groups at the University of Houston site and Hendon | 
_ site. The computer program PILGP1, which uses f-z curves to model pile saapense, 
was used with a shear modulus at small strains to account for the interaction _ 
effects (15). PILGP1 adjusts input single pile t-z curves for the interaction effects 
to account for group behavior. The predicted results were in general agreement 


with observed result 
Conciuoinc Comments 


Any assessment t of the accuracy of the models must ‘consider 
--variability an 
models and our "ability to define the input to ‘the models. Variations in soil 
conditions across a site and in installation details of piles as well as changes © 
in soil properties with time (consolidation, thixotropic and seasonal effects, 
_ and load history) can result i in variations of +30% i in the measured displacement 
at - working loads for ‘supposedly identical piles. For example, load tests on 
four 6-in. (152-mm) diameter timber piles driven 20 ft (6.1 m) into a Bangkok 
clay show a variation of +27% from the mean of the measured pile ire 
at a working load (safety factor of 2) (1). The variation in measured Gees 
capacity can also be large as as shown | by the tests in Bangkok clay | as well - 
_ Accuracy of the measured data is affected | by the care exercised in | taking ‘ll 
a. readings, the resolution and stability of the instrumentation, and the movement — 
of the reference points. emphasis and interest in pile tests 


displacement readings in pile tests may net be sufficient 
to 0 evaluate the > analytic models. The cases s used i in n this study to test the theoretical ic 


- for ‘research: studies where the | accuracy of the measured displacements were ; 


- Predicted and measured load-displacement curves were generally in good 7 
agreement. Improvements are still needed, however, to: (1) Define the distribution 
- and magnitude of the ultimate lo load transfer along the pile shaft and at the 
pile tip; (2) predict the soil movements (and, thus, negative shaft friction) caused 
by consolidation as excess pore pressure, which is induced with pile installation, — 
é dissipates (especially in soft clays); (3) incorporate the effects of soil creep; 7 
(4) evaluate changes in soil poapertios when piles are installed in sands; and 
(5) evaluate group efficiencies, especially if load- -displacement predictions are 
4 The theoretical t-z curves described here can provide reasonable estimates 
>. of the load- -displacement response of single piles and pile groups, oy 
for working loads. These mathematical tools should be used with judgment — 
and knowledge of their limitations. The areas identified as needing wll 
improvements a to ‘these limitations. Our of both 
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__ response of the two piles tested to failure were not in exact agreement with 
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mitations and attributes of the theories will improve as we gain more _ 


- nce with their use 2 and make efforts to compare predicted and measured ay 
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Appenvix 


following symbols are used in this paper: 


= Young’s modulus of soil; 


- at small strains; 


modulus at “small ‘strains and remote from p 


shear modulus at small strains; ry 
= influence coefficient; 62 ines { 


= modulus ratio of soil at pile shaft; Rude 
pie 


Qu = ultimate Q; 
q = stress on pile t tip; 
empirical parameter ofa hyperbolic n 

_ radius of zone of influence; 
undrained shear men th of soi ot 
= shear stress along pile; 


from peak to residual conditions in direct shear 


ratio of soil shear moduli rom = and /; 


shear stress; 


shear stress at pile 7 
o 
= ed max shear resistance; 
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an 
axial pile displacement at pile shaft; tes! 
Az = 
= 


— 
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_ SEEPAGE FROM FREE WATER ABOVE IMPERMEABLE _ 
bus sozmotT to 


Lewis T. Isaacs" and Bruce Hi Hunt? 


on? 


‘Tailings are are the residue from the treatment of ore after te vilken 
have been removed. Coarse, dry tailings are usually stacked on a tailings dump, _ 
ae the fine material is pumped as a slurry to a tailings pond or dam where ~ 
_ the waste material settles. Free water from the slurry, and from mm any runoff — 


from the catchment area a upstream of ‘the tailings dam, may be ponded above a 


the tailings. The dam embankment is usually constructed of readily available 
_ materials which can include spoil from strip mining and dried tailings. Conventional e 
_ practices and typical dams are described by Casagrande and  ooennad (3), Klohn 

(8), and in a number of contributions in Aplin and Argall (1). os 
a Seepage analyses of tailings dams may be done for a number of reasons. 
‘The analyses are used to estimate the amount of water leaving a a tailings structure 

as seepage water, to determine pore pressures for stability analysis, to estimate — 
@ exit point and the seepage velocities. Methods for seepage analysis of tailings | 
Ren are discussed by Osler (9), Galpin (4), Swaisgood —" (13), and | 


‘materials, ‘the conditions for ‘seepage situations in which 
the permeability of the tailings is ; much gremter than that of the embankment, — 


_ the seepage discharge into the embankment 2 and improves its stability. F However, q 
it may not be possible to guarantee this condition, and a designer may have. 

to consider the possibility of the pond water being in contact with the embankment. 

When the tailings are impermeable relative to the bank, the seepage analysis 
_ Tequires an estimate of the seepage discharge from the free water into the — 


iat —'Sr. Lect., Dept. of Civ. Engrg., Univ. of Queensland, Brisbane, Australia. _ 
Reader, Dept. of Civ. Engrg., Univ. of Canterbury, Christchurch, New Zealand. 
Note. —Discussion open until April 1, 1982. To extend the closing date one month, 
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ASCE. Manuscript was submitted for review for possible publication on October 10, 
1980. This paper is part of the Journal of the Geotechnical oe mae Division, Proceedings 


1981. ISSN 0093-6405 / 81 /0011-1563 / $01.00. 
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: dam in such a way that the tailings deposits slope away from the embankment, — 
| 
Le) 
C—O 


- ditch”? solution in which the discharge is assumed to reper half the — 
from a triangular ditch. The discharge for a triangular ditch can be determined 
from the calculations of Vedernikov (5). 
This paper describes a proposed model for the | entry region of the discharge 
from free water above impermeable tailings, and presents an analytical solution — 
based on the model. The results are compared with those of Thomson and 
Rodin (14), and with the ‘half triangular ditch’”’ solution. An experimental 
: investigation is described, and experimental and theoretical results a: are a 


_ When the tailings are impermeable, seepage discharge enters the embankment 
through the part of the embankment face i in contact with the free water. finan 


the embankment 1). In the proposed model, it is assumed that ‘the 

flow in the entry zone is saturated and that this region controls the discharge. 
_ The side boundaries of the entry zone are surfaces where the pore pressure — 
is zero relative to © atmospheric pressure. In the cc core of the flow, the the pore 


ys 


rmeable foundation 


_ pressures sat the top of the entry zone are positive, and transition to unsaturated : 
flow under gravity occurs at the bottom of the entry zone when the pore pressures a a 
become nearly zero throughout the core. In the analytical solution, this boundary — 

a was located at minus infinity, and the results are applicable only when the 
saturated region at the bottom of the embankment is some distance below the 


The problem for which a solution is ‘Sought is shown in Fig. 21 It is assumed aa 


c that saturated, two- dimensional seepage occurs within the region D, CBAD.,, 
that the embankment is homogeneous and i isotropic, and that Darcy’ slaw applies. — 

complex potential function is defined by 
w(z) = o(x,y) + (1) 

: in which z = x + iy; i = V —1; x,y = Cartesian coordinates; ¥ = stream 

function; and the velocity potential, , is defined by 


7 
/ 
| 
i, { 
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2. —Entry- Zone Model in z-Plane 
oe * 


the 
wesc 
gf 


= -k | 


in which k = embankment permeability; p = pressure; p = water mass density; 
and g = gravitational constant. Boundary conditions for and ar are shown 
in Fig. 2, and the corresponding | w- -plane i is shown in Fig. 3. 


ROBLEM SOLUTION AND Resuirs 


‘The problem was by using a is referred to by Aravin 
and Numerov (2) as the “method of boundary-value problems,” 
barinova- ‘Kochina (11) as the ‘ 


‘Thomson & 
ditch’ ‘Solution 


BF 


a 


. 5. —4/kL, ~— and b/L Plotted as Functions of els 
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“this application, the _— “proceeded by m mapping g the w-plane in Fig. 3 into 
the lower half of an auxiliary ¢-plane, which is shown in Fig. 4. Then a 
Riemann-Hilbert problem was formulated and solved to obtain an expression 
for the Zhukovskii function in the t- plane. Thus, the mapping function fo 
w, and the solution for the Zhukovskii function, gave the problem in terms 
of a complex parameter, tins 

The calculated values of the dimensionless « discharge, q/kL, and of the 
dimensionless vertical coordinate of point C, y./L, are plotted as functions — 
of a/m in Fig. 5. It should be noted that the dimensionless Gischarge i is yl 


y anes of the problem. In ‘this approach, the problem is formulated on he 

- complex potential plane (Fig. 3), with x and y as the unknown variables. A 

finite difference method is used to compute numerical solutions. The four solutions _ 
hee are shown in Fig. 5. The agreement between these results, and those 


from the complex variable method, are very good and indicate that the onl 
- solution to the mathematical problem has been obtained. 
The “half triangular ditch”’ solution for g/kL is also shown i in ‘Fig. 5. In 
this solution, C coincides with B and the boundary for | = 0 is a vertical 
_ line. The values of g/kL are less than those calculated by the writers, with» ; 
approximate solution given by Thomson and Rodin is 
which kas! 


—=4sina pred 


- (4) 


sought further details of their work and were informed that 3 was 


q "gives a result about 27% larger than the corresponding value from our «results. 
checks and comparisons Suggest that the correct solution to the 


The symbols defined for ‘the: isotropic : embankment Fig. 2) are “used with 
a bar for the anisotropic embankment. The principal axes of permeability are 
- assumed parallel to the coordinate axes, and k, and k, are the permeability _ 
coefficients for the directions parallel to the x- - and j-axes, respectively. aii ne 


The transformation from the X-f plane to an isotropic x-y plane is given 


q 
7 
| 
velocity has a magnitude equal to the permeability, k, as yo- Plots 
: Of the free surfaces calculated for a/ = 0.23 are shown in Fig.-2, 
| 
r= 


N M R 
OVE 


This transformation leaves vertical dpetetes and p/pg + y the same on | both 

7 a The slope of the upstream face, the length of the wetted surface, and the | 
pe ermeability coefficient for the transformed problem are 7 


az tan” tan a 


<> 


k= (k. k,) 
: ‘The transformation (Eq . 5) and the choice for k (Eq. 8) are such that the _ 
discharge for the anisotropic problem, q, equals the discharge for the isotropic 
problem, g. The value for q (= 4) is found from Fig. 5 for the values of 
and k obtained from Eqs. 6-8. budrom 


An experimental investigation was carried out in the apparatus shown Z 


. elevation in Fig. 6. The dimension normal to this view was 1.09 ft (333 mm). 
The choice of a material to model the embankment was influenced by the condition 

= (1) The capillary fringe was to be as small as possible; (2) the seepage 

o flow had to be laminar so that Darcy’s law was valid; and (3) a ——— 


medium-to-coarse sand which was estimated to be about the upper limit of 
grain size for laminar flow. . The measured capillary rise for this material was 
approximately 1.6 in. (40 mm). The sand was carefully placed in layers of . 
_ approximately 6 in. (150 mm). It was compacted with a vibrator during placing - \ 
; to prevent ‘settlement when the seepage experiments were run. A perspex plate 
was used to model the impermeable layer, and a sealant was applied around 
the edges of the plate to prevent leakage, 
Air trapped in the sand | voids presented considerable difficulties during the 
experimental investigation. The presence of air in the voids near the surface — 
was observed to cause large, approx 50%, increases in the depth for any given © 
discharge because of the reduction in the effective permeability. The results 
presented in Fig. 7 were obtained under conditions in which the air had been 
It was decided that the best values for the pecmsesbiity coefficients of the 
im model embankment would be obtained | from an in situ test. A steady-state test . 
was run in which the ‘upper water level covered all the sand, and the lower 
4 water level | was above the top level of the outlet from the sand. Under these = 


have increased the seepage velocity. The material chosen was a well- =o 


conditions, saturated seepage occurred throughout the sand. The seepage heads 

at 16 positions were measured with piezometers attached to tappings in the 

* wall. The experimental values of seepage head were compared with the 
contours nied seepage head o obtained from three finite element analysis, with k, /k , 
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l 0, 2. 25, and 4.0. The overall between the 
and the analytical results was obtained when the analytical results for k, /k, 


_ = 2.25 were used. For this ratio, actual permeability values were calculated 


= 1.64 x 10° m/s and k, = 0.73 x 10° m/s. These permeability values” 
wail the peoceiiee theory were used to predict the depth- discharge relationship | 
for this embankment under the conditions shown in Fig. 2. The result was 
the line marked theory in Fig. 7. The depth- discharge relationship predicted ; 
from the “half triangular ditch” solution is also shown in Fig. 7. 
__ Experimental observations and measurements were done under the conditions — 
shown in Fig. 6. Dye was injected into the water at the top of the plate, and 
where the surface of the water met the sand. The dye movement confirmed 
the expected pattern of flow in the vicinity of points A and C (Fig. 2), and 
Showed that there was a significant ‘movement of water through t the 


and depth were measured over a range of discharges, and the measured values. 
are shown in Fig. 7. A straight line was fitted to the results for discharges" 
greater than 8 x 10°* L/s, and depths greater than 2 cm. This line was above — 


\ the theoretical line and had a slightly different slope. The fact that the experimental _ a 

' straight line was above the theoretical line and, when extrapolated, did not 
pass through the origin was consistent with the observation of a significant 
‘flowinthe capillary fringes, |} 
While the difference between experimental and theoretical values was relatively 
large at the experimental scale, the amount of discharge through the capillary 
_ fringe would remain relatively constant as the depth increased. Therefore, under = | 
field conditions, where the depth would be approx | m, the difference between ‘ 
the experimental and theoretical discharge values would be about 2% for the | 7 
material used in the model embankment. The difference in slope was probably 4 
caused by errors in the assessment of the permeability values. Thetheory indicates _ q 
that the slope of the discharge versus depth line is proportional to the permeability, a 
and an error of +20% in the determination of permeability values in this experiment | 
would not be unreasonable. Much larger errors could be expected for field | 
situations. Overall, the agreement between the experimental and theoretical results = 
_ was considered to be acceptable, and the writers believe that the proposed 
model and the theoretical results will give satisfactory predictions of the 


most dams, the properties of the materials used i in the embankment 4 


‘there v was as evidence of variation within the ‘model embankment, 
and the evaluation of appropriate permeability coefficients proved difficult. Since 
_ the computed discharge is proportional to the permeability, it may be argued 
that an overly sophisticated analysis cannot be justified in most design situations. 
The theoretical analysis has shown that q/kL lies between 1. 5 and 1.7 over 
most of the range | of a/m. Therefore, for design purposes, the use of g/kL 
is 


| 
} | = the experimental discharge and the analytical discharge and yielded k, : 
q 
ay 
4 
: 
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oe he theory is applicable only if the saturated zone at the bottom of the 
~ embankment is some distance below the level of the free water surface (Fig. 
1). The writers’ opinion is that the | theory may be used when the top of the 
saturated zone is more than about 5d below the free water surface. G%s s 
In general, there will also be some seepage through the railings into the 
embankment If this seepage discharge is judged to be significant, it should 
& estimated by an analysis of the tailings and added to the discharge from 
the free water. The saturated zone at the bottom o of the embankment should - 


be analyzed f for the combined discharge. big at tor Teron orf 


_ Astraightforward application of the Schwartz- Christoffel transformation sual 
the mapping of the the w-plane upon the -plane as 


pig 


q = two- dimensional | flow rate. The real and imaginary parts of Eq. 
9 give expressions for and , respectively, in the lower half of the -plane. 


The Zhukovskii function, = = = 6, i 0) +. 0,(&, is defined 


Since ai ky = = 0 along CD, and AD,, at and : since y cosa — x sina = 0 along 
ABC, the problem for 6(¢) reduces to finding a function that is analytic in— 
‘the lower half of the ¢- t- plane, and that: satisfies the conditions 


8 0) cos a + - 0,(&, 0) sina = (F< 

The function G() is continuous with a discontinuity in 


_ The Reimann-Hilbert problem given by Eqs. 11-12 for 6(t) can b be reduced 
ws a Dirichlet problem by using an auxiliary function defined GD 
. (14) 


¥ in | which a = slope. Thus, Eqs. 
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a(t) F(t) i is in lower t-plane and vanishes as |t| o, and ‘that 


the real axis: 


“the real part of this product satisfies the following boundary conditions along — 


Re F(®)] =0 | <—,1<t 
= \——<&< (15) 


The general | solution of a Dirichlet problem in a mower | 7 mes is ren by 


wr 


FIG. 8.—Integration Contour used to Obtain Eq. 11 


; Polubarino ova-Kochina (11), and the use of Eq. 15 in that result gives sshet — 

The ae Fel denotes the ities: or modulus, of the complex function 
F(t) on the real axis. Eq. 16 gives the solution for 6(t) once the unknown 


constants gand\ have beencalculated. 


& Eq. 16 can be be_rewritten in an alternative way by integrating the function — 
FO In (C + - - oe the contour shown in Fig. 8. This leads: 


i 


SEEPAGE 

— F(t) cota [In(¢+ V —1)+ 


Thus, Eq. 17 be used to rewrite Eq. in the 


Eq. 18 isa alternative to Eq. 16 t— 
Construction of the velocity hodograph for this problem st shows that 


4 e 
abe of Eq. 18 shows that it will ; give the correct result when substituted 


toEq.20if,andonlyif, 


| 


> 
Tole. 


@ = iLke™; and F(t) =| agp —1, This 
“pe 


al 


Thus, Ea. 22 q/ kL to be calculated as as a function of a. 


« Finally, parametric equations of the free surface coordinates along CD, and 
AD. can be calculated by letting ¢ + € for < € < and 1 < < 
“4 oo in Eqs. 9, 10 and 16 after setting w= —ky and w = _ + x — 


— 
| 3 
| 
Eq. 4 
An 
| 
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Evauanion 


One of the most difficult parts of the consisted of 
values for the singular integrals. Infinite limits of integration were converted — 
7 to ‘finite limits LA using the change in in variable € = 1/€’. Then a few of the — 
a Taylor’s series expansion of the 
; integrands. How ever, most t of the integrals were computed by a process known 
as ‘‘subtracting the singularity.’ This relatively little-used method is explained — 
‘Scheid and Kantorovich and and will be shown 


1 1/2 
i+a/a (1 


l+— 
hed a- 0-9” 


d an infinite derivative as § > - 


25, and on the result of integrating the right sides of Eqs. 26-27 can be 
added to obtain the exact — ™ 


GT11 
7 | 
- 
4q 
\FEIGE 
27) 
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FREE WATER 
a 
f@de+ sin a «| @ he = 


arp 
wae 
f and its first both finite at and at € = 1, 
low-order quadrature formula, such as the trapezoidal rule or Simpson’s rule, 
can be used toapproximate the integral in Eq. 28. Furthermore, since the remaining 
“terms in Eq. 28 are the leading terms in a Taylor’s series evaluation of Eq. 
a the approximate value of the integral | in Eq. 28° is only | a small portion — 
of the numerical value of J. In fact, the relative ‘contribution of the integral — 


in Eq. 28 becomes negligible as a — 0, which means that the use of Eq. 28 
leads to an approximate calculation of J that becomes more accurate as the 


of the singularity i in the of (Eq. om increases. 


has been proposed for the entry zone of into an 
from free water above impermeable tailings. A theoretical solution for this model 2 
for two-dimensional, isotropic, Darcy flows has been presented, and the trans- 
formation required if the theory is to be applied to anisotropic embankments — 


— 
q 
| 
| 
| 
: | 
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“are close to the conditions assumed it in the model. = 
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The following sy pmbols are used in this paper: 
3 width of saturated plume; 
depthoffree water; coed 


reaconahly accurate e t rided that i ition — 
H 
‘< 
Casagrande, ., and Mciver, N., esign and Construction of alings ams, | 
— 
— 
= tength of wetted face of embankment, 
p = pressure relative to atmospheric pressure; vo 


WATER SEEPAGE 


_two- -dimensional s seepage d discharge; 


a of velocity vector; 
vertical component of 
complex potextial; 
horizontal distance from origin; a a 
vertical distance from m origin; 


_ Slope of upstream face of embankment, 


= dummy variable defined in 4; 


= Cartesian coordinates; 
— ikz; 


dummy integration variable; ae 
velocity potential, real part of w; and ua ‘ 
stream function, imaginary part of w. | 


- = 
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perscript is used to distinguish values that apply an 
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Discussions 


huabe on any paper presented at any Specialty Conference or other meeting, the Proceedings z 
of which have been published by ASCE. Discussion of a paper/technical note is open to 
anyone who has significant comments or questions regarding the content of the paper / technical a 
note. Discussions are accepted for a period of 4 months following the date of publication 

= of a paper/technical note and they should be sent to the Manager of Technical and Professional 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may — 
be extended by a written request fromadiscusser, 

7 _ The original and three copies of the Discussion should be submitted on 8-1 /2-i in. (220-mm) 


by I1-in. (280-mm) white bond paper, typed double-spaced with wide margins. The length of 
a Discussion is restricted to two Journal pages (about four typewritten double-spaced pages 
of manuscript including figures and tables); the editors will delete matter extraneous to the 
_ subject under discussion. If a Discussion is over two pages long it will be returned for shortening. 
All Discussions will be reviewed by the editors and the Division’s or Council’s Publications 
Committees. In some cases, Discussions will be returned to discussers for rewriting, or they | 
may be encouraged to submit a paper or technical note rather than a Discussion. oe 
Standards for Discussions are the same as those for Proceedings Papers. A Discussion is 
Subject to rejection if it contains matter readily found elsewhere, advocates special interests, 
is carelessly prepared, controverts established fact, is purely speculative, introduces personalities, — : 
. 2 is foreign to the purposes of the Society. All Discussions should be written in the third - 
person, and the discusser should use the term “the writer’? when — to himself. The © 
of the original paper /technical note is referred to as “the author.” A > 
Discussions have a specific format. The title of the original paper/ cochatedl note appears — 
“at the top of the first page with a superscript that corresponds to a footnote indicating the ~ 
_ month, year, author(s), and number of the original paper/technical note. The discusser’s a 
should be indicated below the title (see Discussions herein as an example) with 
The discusser’s title, company affiliation, and business address should appear on the first 
“page of the manuscript, along with the Proceedings paper number of the original paper/technical - 
“note, the date and name of the Journal in which it appeared, and the original author’s name. we 
7" Note that the discusser’s identification footnote should follow consecutively from the original — 
paper / technical note. If the paper/technical note under discussion contained footnote numbers 
q and 2, the first Discussion would begin with footnote 3, and subsequent Discussions would 


= separately to tables and references. In referring to a figure, table, or reference that 
_ appeared in the original paper/technical note use the same number used in the original. : age 


the Publications of ASCE for more detailed information on preparation and submission of 


“ It is suggested that potential discussers request a copy of the ASCE Authors’ Guide at 
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sist benob s vd awode aus Of 
Soi-BENTONITE Siurry Trencu Cutorrs* 


Discussion a by S. A. Jefferis* 


The 2 author has presented some very interesting data on ¢ on the: behavior of of 


soil- bentonite cutoffs. In particular, the concept of filter cake permeability to 
ee 6 ratio, k. /t. provides a most useful design parameter. _However, it 


be used to. and extend the data | 
Using the filtration model presented by Nash , Eq. 3, it can be shown ‘that 
the volume of filtrate, V, obtained up to time, 1, during formation of a filter 


‘in hed C (h, b) i is a ‘eae of h, the formation head and b, a more: general — 
_ variable used to describe the type, condition, and concentration of the bentonite — vi 
slurry. C(h,b) is independent of time and thus can be regarded as a constant 
er tests on a given slurry under a constant formation head. The term v is” 
€ obtained from the theory but is included to represent the small excess — 
flow which sometimes occurs in practice during the first moments of filtration. 
v is often attributed to the flow necessary to establish a filter cake. For times _ 
_ greater than about one minute, Eq. 3 has been found to hold very well in 
ua practice. Differentiating Eq. 3 gives the flow ‘Tate through ¢ the cake bal es 


any time. In this wey, k can be obtained directly from ‘cake 
- formation tests. Thus, it is not necessary to use the full procedure outlined — 

4 by the author of forming the cake for the required time, decanting excess slurry, 

_ and then permeating the cake with water under constant head. The fact that — 
during a cake formation test there is slurry above the cake, and that during 
. permeation there is just water, will not affect the value of k,/t, obtained, — 
as filter cake formation is controlled by the permeation of water through the _ 
 qxisting cake. Strictly, if the author’s procedure is used, the same head should a i 
be used for formation and permeation. If the heads are significantly different . | 

the value of k./t. could be affected. Perhaps the author could confirm that 


_ The validity of Eq. 4 can be checked by replotting the author’s data from 


April, 1980, by David J. D’Appolonia (Proc. Paper 15372), 
Lecturer in Civ. King’ s London, London, England. 


; a a. during a cake formation test, V is recorded as a function of time, C(h,b) ee 
can he fram Eo then he need to cal lt far 


F Fig. 7 as k./t, against Ps ? for constant formation heads. Fig. 16 shows ala 
feplotted data. As far as the author’s actual data points have be 
m™ used. Where data have been obtained from the interpolated curves, the « corre: 
# sponding points in Fig. 16 are shown surrounded by a dotted circle. From 
E Fig. 16 it can be seen that the agreement with Eq. 4 is very good with a: -. 
P slight exception of the 24 hr data. This could imply a breakdown in the theory 
® for long formation times, but this seems unlikely as the fit at 16 hr is stil 
q good and the fit for the one point at 48 hr is also good Se’, eat sear 


While the value of k./t. will obviously depend on the type of bentonite 
= and its concentration, ete., the dependence of k./t, on tf ie may, in part, 
; explain the relatively narrow practical limits for it found by the author (5 x 10° 

to 25 x bd factor for k. /t, a factor « of 

FORMATION HEAD , ft. 


at, 


4 


«FIG. 16. Data Between Filter Cake and 
"Formation Pressure andTime 


(FORMATION TIME , hours) 


25 for formation time (e.g., this could be. changed from 50 hr to 2 ie 
_ The author states that filtrate loss, as measured under the standard American ' 
‘Petroleem Institute conditions, is not a useful design parameter. This is in 
pts with the writer’s experience. Filtrate loss can be related to k./t. 
a paper is currently being prepared which includes this analysis), but even — 
his provides little further information relevant to practice. bak a, 
_ The author’s observations, ‘that the downstream filter cake m may ay be lost when — 
Be the cutoff is put into service, but that the upstream cake is normally stable, — 
are most interesting. During formation (excavation), the cakes may be subjected 4 
j a head comparable to that which they will experience in service (especially — 


: 4. allowance | is s made | for the head loss across the soil- bentonite fill). Also, 


| the bentonite. Thus, from a desk study, one could conclude that the downctreem a 


cake ‘should | be stable. In the cak 
§ be unstable, as any consolidation of the soil- bentonite fill under the applied a 
hydraulic gradient could allow movement and possible rupture of the cake. = 


“1. Nash, Kevin L., “‘Stability of Trenches Filled with Fluids,”’ Journal of the Construction - 
Division, ASC ASCE, Vol. 100, No. CO4, » Proc. Paper 11006, De Dec., 1974, pp. 533-542. . 2 VJ 


Cone PENETRATION IN ‘Som PROFILING 

-f Discussion by! Marius larius Roy* 


the ratio u/ q close to one should be associated with normally coundieael 

clays and lower values occuring at high OCR. Also, they have shown that 

7 u/q,. decreases as the cone angle decreases. Therefore, in order to cxpuineitilly 

_ define a good relationship between u/q. and the OCR, it is appropriate to 
comment on the use of the ratio u/q, in the function “il the geometry uss 
location of the porous element on the pore pressure probe. = Ree te 

‘The piezometer probes developed by Wissa, et al. (31), and 
bear a striking resemblance in their shape (Fig. 9). However, as there is a 
small difference in the shape and in the location of the porous element behind — 
“the lower end of the cone, the writer presumes the magnitude of pore pressures | 
measured with these apparati will be different. On the other hand, if the cone _ 
resistance and pore pressures are not measured simultaneously, it will be more a 
difficult to establish reliably good predictions of the OCR i in 


soil properties of sensitive clays. The basic aim of the piezo-penetrometer is 

_ the simultaneous measurement of cone resistance and pore pressures when the 

apparatus is pushed down at a constant rate of penetration. The instrument 

consists of the well-known Fugro penetrometer, developed in Holland by de 

pain (7), fitted with a porous pore pressure pick up and transducer. The 


has built | an apparatus ‘called ‘piezo- -penetrometer” to investigate the 


point has the classical ‘shape of a cone with an apex of 60° and a base area 


diameter as the base of the conical tip. The pore pressure is measured by 
a pore pressure transducer installed between the cone tip and the load cell. 


7 of 100 mm’. It has a straight cylindrical shaft above the point with the e 
= piezo-penetrometer is equipped with a detachable point in order to a 


“April, 1980, by Mohsen M. Baligh, a Vivatrat, and Charles C. Ladd (Proc. Paper 


Prof. of Civ. Engrg., Université Laval, Québec, Canada. 


= j The authors have presented an interesting approach to describing the clay — 
| 
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the position of pore pressure pick up and to study the pattern of 

pore pressure generation at or above the cone tip. 

_ Five types of tips were used in this study, and details are presented in Fig. el 
i three tips identified as PPA, DPA, and PPS are fitted with a filter 
in the cone tip, while the two other tips, PA or PS and RA are equipped - 
= a cylindrical filter installed immediately above the cone and about one — 


«St Alban test ‘site and reported by Roy, et al. (32), show that maximum n pore 
Pi pressures are measured with the tips PPA, DPA, and PPS as the result of 
3 the strong deformation required to let the piezo-penetrometer penetrate the intact _ 

clay. Pore pressures measured with the tip RA correspond to the new sess 

condition along the shaft after the failure was developed around the conical 
in part of the piezo-penetrometer. Tip PA-PS ‘measures pore seed prevailing, 


he 


on) 


Re. 9. —Schematic Drawing of of Pore Pressure Probe, (a) after Wissa, et al., 


magnitude of the pore pressure ‘to resistance ratio, ulaes was 
calculated by using the pore pressure results obtained with the tips PPA, DPA, 
_and PPS, the tip PA or PS and the tip RA, and the average of the cone resistance. 
Ss (side friction is nearly zero in the St. Alban sensitive clay). The average - 
values u/q, of 0.56, 0.51, and 0. Al: respectively indicate a large v variation in — 
the function of the location of the porous element of the piezo- penetrometer. —_ 
it is well indicated from the results of the authors and the writer that the 
_ magnitude of pore pressures during penetration varies with the shape of the 
_ tip and the location of the porous element on the probe. Therefore, it seems 
important to note that a soil profiling should be established with an appropriate 


tip in order to measure the maximum pore pressures induced during the driving, 


i.e, those measured in the failure zone of the intact clay around the tip. Also, 


at 


FIG. 10.—Details of Piezo-Penetrometer un tant w 


tion test (CPT) presented t to the Tokyo Conference i in 11977. al 


Se. Roy, M., Tremblay, M., Tavenas, F., and LaRochelle, P., ‘‘Induced Pore Ca 
* in Static Penetration Tests in Sensitive Clay,”’ Preprint volume, - 33rd — 
Geotechnical Conference, Sept. 24-26, 1980, Calgary, Canada. 
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Beane OF ON LAYERED SoIs 
na.? M ASt Cc 
se Discussion by A. M. Hanna, M 
— 


: for the numerical estimation of ultimate bearing capacity of footings o on scene 
_ soils. The method incorporates a variety of simplifying assumptions that may 
reduce its validity. The writer wishes to discuss some of these ae. 


1. The case of footings resting on a foundation material consisting of i 
multilayer system has been the object of research over the last two decades. 
The authors’ statement that the bearing capacity of C- soils in a layered system 
has not been investigated seems to be contradicted by the existence of several - 
publications listed in Appendix I. Also, the statement does not bring to the 
attention of interested readers the previous efforts to the stated problem. wy . # 
2. In Assumption 2 the authors define the significant depth as the depth | 
to which the net loading intensity due to structural loading can produce a 
perceptible contribution to settlement or shear stress. The authors proposed — 
Eq. 7 to be used to determine the equivalent significant depth in layered soils. 

pe! the upper layer thickness, z,, was equal to twice the footing width, B, the — 
equivalent significant depth will be equal to ‘ and design of such a foundation 
will depend only on the upper layer soil properties. According to the writer’s 
test results, this could lead to a conservative design if a strong soil layer or 
“below underlies this depth of Z,(= 2B). However, if a weaker layer exists 


= 


below, the bearing capacity could be much lower than predicted by the authors’ 
‘method, especially for the case of strip footings 

§ The writers have some reservations about the use of shear box test results ; 
which have no theoretical verification in predicting the ultimate bearing capacity = 


of square or circular footings. Also, the size of the testing box and footing» 


In conclusion, the writer wishes to suggest that a purely empirical method ‘ 
in geotechnical engineering i is s acceptable provided that the method is 


compared w 

it ber ded fo: r the d 

ecommended for the esigner. 

6 Brown, J. D., and Meyerhof, G. G., ‘‘An Experimental Study of the Ultimate Bearing 


Capacity of Layered Clay Foundations, ” Proceedings, Seventh International Con- 


ference on Soil Mechanics and Foundation Engineering, Mexico, Vol. 2, 1969, pp. 


a “July, 1980, by B. Satyanarayana and R. Garg (Proc. . Paper 15578). 
- > Assoc. Prof., Concordia Univ., — - Engrg. Dept., 1455 De Maisonneuve West, ae 


“Quebec, Canada H3G-IM8. 
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ResorcINoLic Grout | FOR INJECTING SANDY Founpations® 


ot jon by Anth Noel J 


ae resorcinol grout. The writer has several comments to make. — 
of all, the practical history of chemical grouts is rather longer than the authors 
: indicate. Silicate grouts have been in use for 65 yr and the previous use of 
‘resorcinol formaldehyde (RF) ; grouts with both acid and alkaline catalysts was 7 
: described in patents (U. K. Patent 918, 641, and amendment, 1963 to 
Soletanche Ltd.) Alkaline catalysed RF grouts, e.g., MQ 14, from Borden ok) 
Ltd., typically contain 5%-10% of resorcinol, have low viscosities of about 
2 cp at 20°C, and gel times ranging ff6m about 20 min to a few hours. They : 
set to form as very hard, durable, translucent pink solids. Formulation data 
is given in the aforementioned patent specification, but comparisons | of technical < 
_ properties and costs with the acid catalyzed grout described by the authors 
’ “cannot be made because they give no mix compositions. The ome RF grouts = 
certainly meet the aims expressed by the authors, 
q . Acidic RF grouts were used successfully for grouting the Woolwich and Reading 
. beds at Blackwall Tunnel in 1961 (14). The alkaline RF grouts used for pone 
deep rock for mine shaft sinking in North Yorkshire (15) are highly exothermic _ 
7 _and require special batching and proportioning arrangements to prevent premature — 
gelation. Fig. 3 of the authors’ paper suggests that the acidic RF grout is also 
: _ exothermic. Furthermore, Resorcinol is highly dermatitic and requires great care q > 
in handling, as the authors indicate; its unreacted presence | in the ground is 


considerably in increase the ‘strength of deposits with selected chemical 
- grouts and they have been used for this purpose in a number of cases (19). ; 

Nevertheless, designs which seek to use the high strength and durability of © : 
_— grouts such as RF must be considered with great care. Virtually all granular — 


“sols are heterogenous and usually contain zones of low permeability such as _) 
4 fine silts, clayey silts, or clays which cannot be effectively impregnated > lal 
i any form of grout. There have been several basic chemical grout formulations | 
_ previously proposed for use in providing high strength. For example, high 
concentration sodium silicate grouts (16), urea-formaldehyde (U.K. Patent 
876,441, 1960, Soil Mechanics Lt Ltd. ), and oy grouts (U. K. K. Patents 1,019,122 
+ some of th Grout Requirements. .—The authors show that they are aware of hom 
some of the required properties of a satisfactory grout, but they do not clarify © 


restraints which apply in engineering. During the = ‘ 


Binnie & Partners, 


| 
| 


ten number of €f. (17), ‘have placed severe limitations 


n the use of chemical grouts, largely because of their potential for ee) 
groundwater. Two examples may be cited: 


1. Grouts such as AM-9 based on acrylamides have excellent technical — 
_— (5). However, there are well-defined and very low limits for acrylamides = 
‘ in drinking water (21). Many authorities will n no longer accept the risk of pollution 7 = 


by si such grouts and their use has declined greatly i in recent years. icant Gee... 

2. Lignochromes first used for engineering purposes in 1957 (19) are are prepared 

_ from lignosulphonates, sodium dichromate, and inorganic catalysts. They, too, 

are excellent and cheap chemical grouts and were used successfully for several é 

major projects (13). Their use has since been in part discontinued because of = 

fears of pollution by chromium. The acceptability of RF grouts must ; 
some doubt for similar reasons. 


doubts about their durability. A further class is based on tannins and 
which “commercial examples ai are manufactured in the form of “Geoseal” 
and “‘Terranier.’’ These grouts tend to contain small particles which limit their 
use in fine-grained soils unless costly clarification procedures are used. A third 

class is the RF grouts which are also costly. 
‘aa Mechanisms of Ground Strengthening.—The reaction of resorcinol with formal- _ 
dehyde under alkaline - conditions | peoduces continuous hydro gels w which fill 
the pore space of granular soils. The resultant product may be regarded as 
= a matrix of set grout densely filled with sand. Indeed, the same description — 


can be used for the majority of chemical grouts because they also set to form : 
_ There are two types of gi grout known to the writer which confer strength — 
on the ground by depositing adhesive material at the contacts of the soil particles: 4 
ny They are (1) the more dilute urea-formaldehyde (UF) grouts; and (2) the special — 
nonaqueous epoxy grouts of low viscosity described in the patent literature — 
of oil well drilling. From the authors’ description, it would seem that the acidic = 


Further Comments.— —No hydro gel chemical grout can endure 


ws inadvisable to use RF ‘grouts in 1 these conditions. The “hydrochloric acid RF j 
; — described by the authors has a pH of 1.5-2.5. The writer would expect 


13. Black, J. J. -Kelland, J.D. , ““Cominco’s Saskatchewan Potash Shafts,” 


of the Ninth Commonwealth and Mining Congress London, England, Vol. 1, Mining _ 

14. CaronC., Deslisle J. P., et al., ‘ ‘*Resin Grouting with Special Reference to the Treatment 

wi Woolwich and Reading Beds at the New Blackwall Tunnel , & Drilling a 
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15. ‘Cleasley, Pearse, G. E, et a. “Shaft at Boulby Mine 
Potash Ltd.,”’ Institute Mining & Metallurgy, Vol. 85, 1975, pp. A7-A28. seu orl no 
. Clough, Wayne G. ., Kiick, W. M., et al., “‘Silicate-Stabilised Sands,’’ Journal of 
the Geotechnical Engineering Division, ASCE, Vol. 106, No. GT1, Proc. Paper <i 
. European Standards for Drinking Water, Second Edition WHO, 1970. 
3. Golden, J., Ouellete, R. P., et al., ‘Environmental Impact | Data Book,” Anne Arbour 
19. 19. James, A. N., and Neelands, R. J., ‘‘Formulation and Selection of Chemical Grouts 
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The discusser has 
alkaline RF 


grouts possess 5 low | viscosities of ‘about. . cp at 20° C and py times ranging, 
- from about 20 min to a few hours, they exhibit low gel strength and 5%-10% 
_ synerisis, which is not with acid RF grouts developed by the writers. Moreover, 


a characteristics of acid catalyzed RF grout are 30%-40% less a a 


AS the writer’s work is limited to academic interest, the efficacy of grouts 


_ for strengthening granular soils with zones of low permeability, such as Lead 


silts and clayey silts in the field, has not been examined. 
_ The effect of leached water from the grouted sand ai with resorcinolic 


a grouts on fish life and the potability of ' water could be examined to conform > 


_ While the corrosion effect on the conventional grout pumps should be examined, 


currently, anticorrosion chemical pumps are employed in 
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published by the American National Standards Institute and to the Authors’ Guide to the Publications ; 
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_ 11. A summary of approximately 40 words must accompany the paper. 
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